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SUBJECT Program for the Solution of Multipoint 

Boundary Value Problems of Quasilinear Differential Equations 


The program is designed to provide numerical solution of multipoint 
boundary value problems governed by a system of up to 20 first-order 
ordinary differential equations. This system of equations will be 
subject to a set of up to 20 boundary conditions. Linear equations 
are solved by a method of superposition of solutions of a sequence 
of initial value problems. For nonlinear equations and/or boundary 
conditions , the solution is iterative and in each iteration a problem 
like the linear case is solved. A simple Taylor series expansion is 
used for the linearization of both nonlinear equations and nonlinear 
boundary conditions. The perturbation method of solution is used in 
preference to quasilinearization because of programming ease, smaller 
storage requirements, and experiments indicate that the desired con- 
vergence properties exist although no proof of convergence is given. 
The use of this program for nonlinear differential equations is also 
applicable to the inverse problem of; given a model and an output, 
determine the system that gave this output. These are often called 
"parameter estimation" or "systems identification" problems . ' 
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PROGRAM DESCRIPTION GUIDE 


IDENTIFICATION 
Program Name - MINI QUASI 

Programmer's Name - V. J. Marks and B. G. Hall 
Date of Issue - 4/9/73 


GENERAL DESCRIPTION 


The program is designed to provide numerical solution of multipoint 
boundary value problems governed by a system of up to 20 first-order 
ordinary differential equations. This system of equations will be 
subject to a set of up to 20 boundary conditions . Linear equations 
are solved by a method of superposition of solutions of a sequence 
of initial value problems. For nonlinear equations and/or boundary 
conditions, the solution is iterative and in each iteration a problem 
like the linear case is solved. A simple Taylor series expansion is 
used for the linearization of both nonlinear equations and nonlinear 
boundary conditions. The perturbation method of solution is used in 
preference to quasilinearization because of programming ease, smaller 
storage requirements, and experiments indicate that the desired con- 
vergence properties exist although no proof of convergence is given. 
The use of this program for nonlinear differential equations is also 
applicable to the inverse problem of; given a model and an output, 
determine the system that gave this output. These are often called 
"parameter estimation" or "systems identification" problems. 

USAGE AND RESTRICTIONS 


Machine and Compiler Required - Xerox Sigma 3 

- ANSI FORTRAN 

Peripheral Equipment Required - Card Reader and Line Printer 

Approximate Amount of Memory Required - 13B3 liexidecimal (5043 deci- 
mal) 


PARTICULAR DESCRIPTION 


Equations Used and Derivations - See Reference 1 

Definition of Terms Used - See Reference 1, Appendix A and B 

Detailed Description - The program is designed such that 

a broad class of problems can be solved with a minimum amount of pro- 
gramming effort by the user. The on3y programming required by the 



user is that which is involved with the description of the governing 
equations for a specific problem, and this task is completed within 
the framework of a subroutine furnished with the program package. 

All other pertinent information is communicated .to the program by 
means of input data cards. 

After setting up COMMON and DATA blocks, the main program begins exe- 
cution by reading data to obtain values for certain key control para- 
meters, boundary condition information, and estimates of initial values 
of the solution. This data is then reduced as follows: the boundary 

conditions are arranged in order of their occurrence with respect to 
the independent variable; the limits on the independent variable are 
changed if need be to span the boundary conditions; any initial con- 
ditions which were input as boundary conditions are placed in the 
initial value vector where they supersede any previous values , For 
each exact initial condition, the order of the system to be solved is 
reduced by one until the required number of solutions to be used in the 
superposition scheme has been determined. The perturbation factors are 
then checked to ensure that none of the perturbations are null, and the 
required matrix dimensions are calculated. Finally, an echo check of 
initial conditions, input data, and key program parameters is accom- 
plished by outputting appropriate values . 

Iterations 


Next, up to a certain number of iterations are taken in an attempt to 
solve the problem as input. ' The initial values of the best estimate 
of the solution at each iteration are used to generate both the "prev- 
ious approximation" and the particular solutions to be used in forming 
the next approximation. These solutions are started from initial values 
which are perturbed' from those of the previous approximation in the 
manner already discussed. The integration of these initial value prob- 
lems is undertaken by a fourth order Runge Kutta scheme especially adapted 
for the data structuring employed in QUASI. Integration starts at the 
minimum value of the independent variable and stops at the last prescribed 
boundary condition during intermediate iterations and at the maximum value 
during the final iteration. Output of integration results is controlled 
by the user through an input parameter. 

The superposition of these solutions is implied by interrupting the inte- 
gration whenever a boundary condition is encountered and setting one of 
the algebraic equations into a storage area for future use in determin- 
ing the superposition constants. When all boundary conditions have been 
so set,’ a linear algebraic system is solved to obtain the solution as 
desired (exact fit or best fit in least squares sense). Following the 
determination pf the constants, the best estimate of the initial values 
is .updated according to the desired method and if convergence has not 
occurred, within the limiting number of iterations, a new iteration is 
begun. Convergence occurs whenever all of the constants of the perturbed 
solutions are smaller in absolute va] ue than a preassigned value . 
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e'. description of input 

(Extracted from Ref.l) 

In order to utilize the program the user must successfully communicate a 
variety of descriptive- and control information. There are two basic chores 
involved. The first consists of programming the system of equations into 
the structure -of one of the FORTRAN subroutines . The second task involves 
the input of program control parameters, boundary value information; and 
other data necessary--or desirable to a successful program execution. 

Nonlinear System Programming of Equations 

The way in which the nonlinear system may be described in the RHS sub- 
routine is as follows 

Denoting the previous approximation to the solution of the nonlinear system 
y = g(y,t) 


as 


V 


the linearized system can be written 
. dg(y ,t) 

7 ‘ g( V 10 + ly ' y ° ) 


dg(y ,t) 

where denotes the Jacobian of the nonlinear system 

oy 

■evaluated using the approximation y . Evaluating the indicated diff- 
erences in current and previous approximations leads to a system of 
equations 




>t) 


V a8 i ^ y o ,t) |> ft) - y (t)] 

U V ° °j 

J 


These equations can be programmed as above using the vector Y for 
the current solution and the vector YO (or Y0) for the previous approx- 
imation. .Of course, the values of g anddg../dy_ are to be computed 
using the previ ous approximation. * 1 i 

Linear System 

The actual linear equations are placed in the RHS subroutine. 

No further statements are required. ' 

1. Data Input 

The input data is read by the program in four groups: integer control para- 
meters, real control parameters, boundary value information, and initial 
conditions and related data. The input for each section is terminated by 
a blank record and the procedure is so structured to allow an arbitrary 
ordering of records in the input list and an arbitrary amount of data to be 
read in any one of the four groups. This input pi'ocedure is particularly 
convenient for multiple data sets for the same system of equations, because 
only that information to be changed need be included in the input list. 



k 


The integer and real control parameters have been included in two 
COMMON blocks, and input of values for these parameters is done by- 
reference to the order of occurrence of a particular element in a 
COMMON block. This practice greatly simplifies input commands. 


Integer Control Parameters 


The elements of the integer COMMON block have been EQUIVALENCE! to 
elements of a vector of integer variables occupying the same storage 
location and can be referenced by name or by lineal subscript location 
within-..the vector. These elements are input by entering the vector 
subscript corresponding to the element location and the desired value 
of that element in the first two fields of a record. Listed below are 
the most commonly used 'input integer variables with the vector subscript 
as the key variable. The variables are all given a default value. If 
the value is nd£ listed, it is zero. In most cases, the zero value is 
an ideal value . The - format for integer input is 


FORMAT (215) 


ELEMENT 

NUMBER DESCRIPTION 

1 output control parameter. Zero default. 

= 0 gives minimal output. This includes the initial conditions of 
each iteration and the superposition constants . On the final iter- 
ation the unperturbed or nonlinear solution is output. 


= 1 adds to the above output the initial condition change that would 
be taken if the change is unconstrained, the cartesian length or 
norm of this vector and the. dot product of .successive unit" vectors 
parallel to the changes. The matrix in which the boundary conditions 
and results of the boundary condition operators on the perturbed 
solutions is also output . 

= 2 adds to all the above output, the pivotal elements used in 
solving for the superposition. 

= 3 adds the unperturbed or nonlinear solution to the above output. 

= 4 adds the perturbed solutions to the above output. 

2 =0 for nonlinear problem. Zero default. 

= 1 for linear problem. 

3 See Appendix D for explanation. Zero default. 

h = maximum number of recurrence iterations . 

Default value = 5* 
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ELEMENT 

NUMBER 

5 

6 

7 

8 
9 

10 

11 


DESCRIPTION 

= 0 keeps a constant maximum value for the norm of the change of 
the initial condition vector. Zero default. 

0 allows the maximum value for the norm of the change of the initial 
condition vector to be reduced to the smallest norm' encountered. 

= total number of boundary conditions to be imposed on the system. 

Must he input except on initial value problems . 

= number of dependent variables in the problem which have non -null 
derivative expressions. Must be input because default is meaning- 
less . 

= number of dependent variables which do have null derivative expres- 
sions; i.e., constants to he identified. These are the last elements 
of the state variable vector. Zero default. 

output control parameter which gives frequency at which output of 
information about the solution is to occur. It refers to the nominal 
number of integration steps between regular output points . It should 
be noted that output will also occur at points where boundary con- 
ditions have beep specified. Default value =5. 

5 L 0 causes the original initial conditions to be saved for use with 
subsequent data sets. 

= 0 the subsequent data sets will use the last information [then avail- 
able for initial conditions. Zero default. 

used to signal computation and output of the Wronskian of the solution. 
Zero default. 

< 0 no computation 

= 1 computation only at end of integration on final solution. 

> 1 computation each time a boundary value is encountered during 
solution. 


The remaining elements in the integer COMMON block are established during execu- 
tion and heed not be input . A blank record should follow at this point to term- 
inate integer input. 



6 


Real Control Paratneters- 

manner- as those 
field specifi- 
is listed below 
real parameter 

FORMAT (15,G15.7) 


The parameters in the real COMMON block are input in the same 
in the integer COMMON block (except, of course, for different 
cations). The vector subscript is again the key variable and 
with the corresponding parameter description. The format for 
input is 


ELEMENT 
NUMBER ' 

1 

2 

3 


h. 


5 

' 6 


7 , 8 
9,10 
ll' 

12 


DESCRIPTION 


the nominal integration step size to be used' — default value - 0.1.' 

a standard perturbation factor to be applied to all unknown initial 
value estimates unless otherwise specified — default value = 1.2. 

the maximum value of the "length" of the change in the initial value 
vector- -if input as zero, will override 5th element of integer 
COMMON block being non-zero and allow unrestricted changes. Zero 
default . 

the initial value of the range of the independent variable . It is 
set to the smaller of (zero or the input value) or the smallest time 
of a boundary condition. 

the final value of the range of the' independent variable . It is ; 
set to the larger of (zero or the input value) ’or the largest time 
of a boundary condition. Must be input for initial value problems. 

a scale factor for the determinant of the coefficient matrix used 
'in calculating the superposition constants . Default value is unity. • 

small parameters which need not be input. 

parameters which should not be input. 

a small parameter used to check on ill conditioning while solving 
for 1 superposition constants. Default l.E-20. 

convergence test parameter. Default l.E-8. 


13 parameter used in checking for nearly null initial values . 

Default l.E-3* 


Boundary Value Information 

For each boundary condition imposed, a record consisting op five fields is 
read. The first two fields .are integer values, the next two fie.l da 'arc real 
values, and the last field is another integer value; The format used is 


’FORMAT (215,2015. 7, is) 



FIELD 


DEFINITION 


1 an identifying- integer for the boundary condition — these values 
should range from one through the' number of boundary conditn ons 
used, with no repetition--they will -later be arranged in order 
according to the ascending values of the independent variable at 
which the boundary condition occurs . 

2 the element of the state variable vector on which the condition 
is imposed. See Appendix D for general linear and nonlinear 
boundary conditions. 

3 the value of the independent variable at which the condition is 
imposed. 

4 the numerical value of the boundary value . 

5 a flag indicating whether the boundary condition is to be met in 
an exact sense or a best-fit sense-- 

= 0 implies exact . 

>0 implies least -square. 

After including one ' such record for each boundary condition in this group of 
data, a blank record should again be used to terminate input. 

Initial Conditions and Related Data 


QUASI requires that initial conditions, either exact or estimated, for each 
element of the state variable vector be available at tfte beginning of a solu- 
tion. It is often the case that some boundary conditions are specified at an 
initial value of -the independent variable, in which case they automatically 
are placed in the initial condition vector. For each element of the state 
variable' vector for which this does not occur, an. estimate of the initial con- 
dition should be input at this time. See input of integer parameters, 10th 
element and SKIP, Appendix B of Reference 1. Each initial condition to be 
input is included on one' record consisting of two integer fields followed by 
four real fields read according to: 

FORMAT (15,G15.7, 3G10.3) 

FIELD ' DEFINITION 

1 integer indicating the element of the state variable vector. 

2 - 0 implies an unknown initial condition without bounds (see 
fields G) 

- -i 1 implies <m unknown iuilial condition with hound;:. 

-.1 i in | t I i«\; an exacl, ini I. in I .miidi l.imu If i.lvis i:. l-lio ■vu-a;, 

(.ho fields, J l , and <> an> ignored. 


ORIGINAL PAGE IS 
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.FIELD DEFINITION" (cont'd) 


3 


the floating-point value of the initial condition. 
Default value = 1.0. 


1 the value of an associated perturbation factor, if different 

from standard value or zero. 

5', 6 if one of these parameters is non-zero, the larger of the two 

values will be used as an upper bound for the initial condition 
and the other value will be a lower bound and if the second field 
is null, it is then changed to +1. 


A blank record is once again used to terminate input of this group of data. 
This concludes input of a complete data set. If subsequent- data sets are 
incomplete, i.e., with one or more of the four data groups empty, simply 
include a blank record as a terminator for. the empty set and continue as 
before i • ‘ 


2. Control and Program- Cards - (begin Col. l) 

!J0B QUASI • 

IF0RTRAN 

(Source Cards - See Appendix A) 

IE0D 

I0L0AD 

1$R00T 256, ,G0 
:$mp 

IEND 

1XEQ 

‘(input Data) 

:e0d 

F. DESCRIPTION OF OUTPUT 

Printer Output - The iteration number and solution are given for each 
boundary value. See Appendix A for sample output. 

G. INTERNAL CHECKS AND . EXITS 


Program exits when input conditions for convergence arc mot. 

H. INDEPENDENT SUBROUTINES 


None 



I. 


SYSTEM SUBROUTINES 


y 


No special subroutines . 

J. COMPLETION OR FINAL CHECKOUT DATE - U/9/73 

Source listing and sample cases are given in Appendix A. The first 
sample problem is a duplication o P the first sample problem discussed 
in Reference’!. In this problem the boundary conditions are nonlinear 
and are fitted in a least square sense. The .programming of the non- 
linear boundary condition equation is done in subroutine' NLBC . 'The 
second example problem solves the same system of equations appearing 
in subroutine RHS, but linear boundary conditions are used and the 
subroutine NXBC is ignored'. In the second example problem, the linear 
boundary conditions are also fitted in the least square sense. 
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APPENDIX a 


SOURCE LISTING 
AND 


SAMPLE CASES 



_L C. 


X 

2 

3 

4 

5 ■ 

6 
7 
B 
9 

19 

it 

1.2 

13 

14 

15 

16 

17 

18 
19 
29 
21 
22 

23 

24 

25 

26 
2 ? 
28 
29 
39 

31 

32 

33 
3$ 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 
47' 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 
58. 


Q,** ************************************************************ ********* 

q ****.*******^ ****** *********** QUASI ***************************** 

Q** ******** ******************* ***************************** 

c *********************************************************************** 

c • 

INTEGER Q8V.#XACTIC#SKIP# IBIG( 135)#XACTBV 

REAL RBIG(204 )#LOWER#SVEC( 462 ) # LF# IC#IIC#C( 462 ) # Y < 20# 21 ) 

COMMON IIC(20)#1C< 20>#YVEC( 420)#YNL( 80 ) # LF < 20 > # S C 21# 22 > # 

1 DELTj SPTRB# PNORM# - TSTART# T STOP# 

1 DET# SMALL# ZERO# DETT# TRACE# 

2 DZER0# C0NV# DPTRB# DTEMP# DN1 # 

3H#RN1#RN2#T# 

4TAVG#TL#TN# TBMT#TEM#TEMP# TTEM#TBV( 21)#BV( 20)# 

5 UPPER ( 20 ) # LOWER { 20 ) # PTRB ( 20 )# AUX ( 10 ) # CQBV! 20 )# 

5 Dl C I 20 )# YPTRBT 20 ) # DELTP# DELTM# TNC1RM# S N 8 R M / AN0RM#TBVK#D0T# 

1 I6UT# UN# IQBVS# LITER# N0RMRD# 

1 NB V# ■ NEQ# N C N # NWRITE# SKIP# 

l' • ITRACE# LIMIT# IVARY# IAF(2)# 

1NVECT# IXTRA# 

1 JXTRA#KXTRA# N0UT# NL#NSPACE# NRCD# 

‘ 2NRSD# NR YD# I# 1 1# IC0UNT# I Nl# I N2# INTKEY# IPUNT# IQ# IS# ISUB# I TEMP# I TER# 

3 J# JJ# JQ# JSUB#K#KK# K0UNT# K0UT#KQ# KTEM#MITER#NBVM# NC0L#NRANKC# 

4 NR0WS#QBV(2O)#XACTBV(21)#XACTIC(2O># IQBV (20)#NR#IG0#NE#NI 
EQUIVALENCE ( YVEC # Y ) # f S # S VEC ] # ( YVEC#C ) # „( I OUT# IBIG)# ( DELT# RB IG ) 

C****WARNING# CHECK COMPILER F0R THIS OPTION**** 

C ****InPUT INTEGER PARAMETERS**** 

C****WARNING# CHECK COMPILER FOR THIS 0PTI0N**** 

WRITE- (.102# 15 ) 

WRITE ( 6# 1 5 ) 

15 FORMAT (/' MINI “QUASI MODEL'/ 

* ' REFER TO GE-AGS USER GUIDE TIR 74 1 4MED-30XX » // ) 


CALL INIT 
1000 CONTINUE 

CALL READIT 
CALL S0RTBV 
CALL NTLIZE 
CALL ECH0IN 
5 CONTINUE- 
CALL ESTABL 
CALL NGRATE 
I VARY=*4 
7 CONTINUE 
L I M I T= s 4 
' CALL ONETEN 

GO TO ( 17# 18# 19 ) # IVARY 

17 CONTINUE 
CALL WRITES 

20 GO TO < 7# 18)# I GO 

18 CONTINUE 
CALL STORES. 

IF IK-NBV ) 7#7# 19 

19 CONTINUE 

IF ( ITER-MITER ) 25# 1000# 1000 
25 CONTINUE 
LIMIT-1 
ll CONTINUE 
CALL WRJTES 

GO TO ( 12#5# jlOOO)#LIMlT 




. 1.3 


59 

12 

CONTINUE 

60 


OETT = DET 

61- 


CALL GJRWLS 

62 


LIMIT = 2 

53 


IF ( IPUNT-2 > 13; 11;13 

64 

13 

CONTINUE 

65 


CALL NEW I C 

66 


IF (J) 14;11;11 

67 

14 

CONTINUE 

66 


CALL MODIFY 

69 


GO T 0 13 

70 


-END 


71 SUBROUTINE INIT 

72 C 


73 
7 4 

75 

76 

77 

78 

79 
89 
81 
82 
83 
8# 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

96 

99 

97 
9B 
99 

100 

101 

102 

103 

104 

105 
10 6 
10-7 
108 

109 

110 
111 
112 

113 

114 

Ml 

116 

117 


INTEGER 0BV;XACTIC;SKIP;IBIG<135);XACTBV 

REAL RB IG ( 204 ) * LOWER/ SV£C ( 4 6 2 ) » L F; I C > I I C ; C < 4 62 ) ; Y { 2 0/ 2 i ) 

COMMON IIC(20);1C<20);YVEC( 420);YNL< 80 ) ; L? (20);8(21;22); 

1 DELT; SPJRB; PN0RM; TSTART; TST0P; 

1 DET, SMALL; ZERO; DETT; TRACE; 

2 DZERO; CONV; DPTRB; DTEMP; ONI; 

3H; RN1 ; RN2; T ; 

4TAVG; TL;'TN; TBMT; TEM; TEMP; TTEM;TBV< 21>;BV< 20); 

5 UPPER( 20 ); LOWER! 20 ) ; PTRR (20);AUX<10);CQBV<20); 

5 DIC{20 ); YPTRBt 20 );DELTP;OELTM;TNQRM;SN0RM;AN0RM;T8VK;D0T; 

1 IOUT; LIN; IOBVS; LITER; N0RMRD; 

1 NBV; NEQ; NCN; NWRITE; SKIP; 

1 ITRACE; LIMIT; I VAR Y; IAF(2); 

1NVECT; I XT R A ; 

1 JXTRA;KXTRA;N0UT;NL;NSPACE;NRCO# 

2NRSD; NR YD ; I ; 1 1 ; IC0UNT ; INI ; I N2; I NTKE Y; I PUNT; 10; IS; I SUB; I TEMP; I TER; 

3 J; JJ; JQ; JSUB; K; KK; K0UNT; K0UT;KR; KT EM ; M I T ER; NB VM; NC0L; NRANKC; 

4 NR0WS;QBV(20);XACTBV(21 );XACTIC(20); ieBV(20);NR; IGO;NE;Nl 
EQUIVALENCE ( YV£C;Y ) ; ( S;SVEC) ; ( YVEC;C ); ( I OUT; IBIG);<D£LT/RBI6) 

DATA I DUM; I DUN/ < 

DO 10 M = 1 ; 1 35 
10 J8IG(H)*0 

RBIGT 1 ) =1 .E«l 
DO 20 I *=2;204 
20 RB I G ( I ) =0 » 

00 30 I e= 1;20 
30 I IC< I ) = • 1E1 
NRSD *» 21 
NRCD * 21 
NRYD - 20 
L I TER =5 
DET=. 1E1 
SPTRB® • 12E i 
SMALL** IE-9 
ZER0=»5.E-6 
DZER0~ 1 * E *> 2 0 
DPTRB = l#E-3 
CONV * l.E-8 
NWRITE«5 
I AF ( 1 ) * I D U M 
IAF(2 >=IDUN 
RETURN 

END 

SUBROUTINE READIT 


C 



118 INTEGER QBV/XACTIC/SKIP/IBIGI 1351/XACTRV 

11-9 REAL RBIG(2O4)/L0WER/SVEC( 46£ ) /LF/ I C/ 1 1 C/ C ( 462 ] / V ( 20/ El ) 


120 

COMMON 

IIC(20)/IC(20)/ 

YVEC ( 420 ) / YNL ( 80)/LF(20)/ 

S ( 21/ 22 > 

121 

1 

DELT/ 

SPTRB/ 

PNORM/ 

TST ART/ 

TST OP/ 

122 

1 

DET/ 

SMALL/ 

ZERO/ 

DETT/ 

TRACE/ 

123 

2 

DZER0/ 

CONV/ 

DPTRB/ 

DTEMP/ 

DN1 / 

124 

3H/RN1/RN2/T/ 






125 

4TAVG/TL/TN/TBMT/TEM/TEMP/TTEM/TBV( 21 ) 

/ B V ( 20)/ 



126 

5 UPPER( 20 )/ LOWER ( 20 )/PTRB (20 ) 

/ AUX ( 10 )/ 

CQB V ( 20 ) / 



127 

5 DIC( 20 )/ YPTRB ( 20 ) / DELT P/ DELTM/ TNORM/ SN8RM/ ANORM/ TBVK/DOT/ 

128 

1 ^ 

I OUT/ 

LIN/ 

IGBVS/ 

LITER/ 

N0RMRD/ 

129 

1 

NBV/ 

NEQ/ 

NCN/ 

NWRITE/ 

SKIP/ 

130 

1 

ITRACE/ 

LIMIT/ 

IV ARY/ 

I AF ( 2 ) / 

131 

1NVECT/ I XT R A / 







132 1JXTRA/ KXTRA/N0UT/ nl/ nspacE/ nrcd/ 

133 2NRSD/ NR YD/ 1/ I I /.I COUNT/ JN1 / IN2/ I NT KEY/ I PUNT/ 10# IS/ I SUB/ I TEMP/ ITER/ 

134 3 J/ JJ/JQ/ JSUB/K/KK/ K0UNT; KOUT/KQ/ KTEM/ MITER/ NBVM/NC0L/NRANKC/ 

135 4 NR0WS/ QBV( 20 )/ XACTBV ( 21 > /XACTICI 20 )/ IGBV < 201/NR/ IGO/NE/NI 

136 EQUIVALENCE < YVEC/ Y ) / ( S/ S VEC ] / ( YVEC/ C > / < I OUT/ 1816)/ (DELT/RBIG) 

137 WRITE (6/312) 


138 

1 

REAO ( 5/ 300 ) INI/ IN2 

139 


I F ( INI ) 230/ 5/ 7 

J40 

5 

WR I TE ( 6/ 302 ) 

1 4 t 


GO TO 11 

142 

7 

WRITE (6/301 ) INI/ 1N2 

1<*3 


I B I G < INI ) *IN2 

144 


GO TO 1 

145 

C****INPUT REAL PARAMETERS**** 

146 

11 

READ(5/305> IN1/RN1 

147 


I F ( INI ! 15/15/17 

148 

15 

WR I TE ( 6/302 ) 

149 


GO TO 20 

150 

17 

WRITE(6/306)IN1/RN1 

151 


IM100=IN1-100 

152 


I F ( IM100 ) 18/18/ 19 

153 

18 

RBI 6( INI ) =RN1 

154 


GO TO 11 

155 

19 

AUX! IM100)=RN1 

156 


GO TO 11 

157 

C** ** INPUT BOUNDARY VALUE INFORMATION**** 

158 

20 

IQBVM = I QB VS 

159 

21 

READ ( 5/ 300 ) 1/ IN1/RN1/RN?/ IN2 

160 


I F ( 1)901/901/905 

161 

901 

WRITE! 6/302 ) 

162 


GO TO 31 

163 

905 

WRITE (6/301 ) 1/ IN1/RN1/RN?/ IN2 

164 


IF ( INI ) 906/ 907/ 907 

165 

906 

lF(lNl+100)907/907/908 

166 

907 

QB.V (I ) = IN 1 

167 

908 

TBV ( I )=RN1 

168 


B V ( I ) =RN2 

169 


XACTBVI I ) =1N2 • 

170 


I F < INI ) 909/ 21/21 

171 

909 

IF(IN1+100). 9 10/9J0/22 

172 

910 

I F ( QB V ( I ) >911/22/ 22 

173 

911 

I N 2 = Q B V ( I )/100 

174 


' I F ( IN1"(0BV(I )~IN2*100) )?2/912/91? 

175 

912 

I QB VS=" IN2 

176 

22 

IF ( IN 1+ 100 ) 851,851/850 
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8SG QBV ( X > = «»{ < IQBVS + 1 >*100-IN1 ) 

851 INI * -Q0V{ I )+(QBV( n/ioouioo 
IF< INI ) 913, 21 #913 

913 DO 23 1*1, INI 
IQBVS « 1 GB VS + 1 

RE AO t 5, 305 ) IQBV < IQBVS 1 , CQBV ( IQBVS ) 

I F I IQBV ( IQBVS >> 23, 914, 23 

914 I QB V { IQBVS ) =1 

23 WRITE { 6,306) I QB V ( I QU VS ) , CQB V ( I QBVS ) 

IF I IQBV S« IQBV M> 24, 21, 25 

24 I QBVS* I QB VM 
GO TB £1 

25 I QB VM* T QBVS 
G0 TB 21 

C****INPUI INITIAL C0NDITI3NS AND RELATED DATA**** 

31 READ ( 5, 3 1 D ) I N 1 , IN2, DN 1 , RN 1 , RN2, RN 3 
I F ( INI 132,32,915 

32 WRI TE ( 6, 30£ ) 

G0 T0 34 

915 WRITE { 6 , 311 H N1 , I N2, DN1 , RN1 , RN2, RN3 

1 1 C 1 INI ) « ONI 
PTRB {INI) « RN1 
XACTIC (INI) = IN2 
IF< A8S{ RN2)-DZERfl >916,916,917 
91b IF(ABS(RN3)»DZER0 >31,917, 917 

917 IF(RN2-RN3) 918, 918,33 

918 T EMP a RN2 

RN3 « TEMP 

33 UPPER < INI ) s RN2 
LBWERt INI > * RN3 
XACTIC (INI) » 1 

G8 TB 31 

34 RETURN 
230 ST8P 

300 F0RMAT (2I5,2E15» 7, 15 ) 

301 FORMAT (16H INPUT DATA CARD, 2 1 5, 2E1 5 - ?i 1 5 ) 

302 FORMAT (15X,16H BLANK DATA CARD, // ) 

305, FORMAT { 15, El5»7) 

306 FORMAT (16H INPUT DATA CARD, IS, E15,?J 

310 FORMAT ( 2 1 5, E 15 , 7 , 3E 1 0 . 3 ) 

311 FORMAT (16H INPUT DATA CARD, 215, 4E1S • 7 ) 


FORMAT ( 1 H 1 ) 


SUBROUTINE S0RTBV 


C****SORT BOUNDARY VALUES**** 

'INTEGER QBV, XACTIC, SKIP, IBIGl 
REAL RBI.G (204 ), LOWER, 

COMMON IIC{20),IC(20>, 

l DELT, 

1 DET, 

2 DZER0, 

3H, RN1, RN£, T, 

4TAVG,TL,TN,TBMT,TEM,TEMP, TT EM 
5 UPPER ( 20 ), LOWER ( 20 ),PRTB( 20) 
5 DIC( 20 ),YPTRB( 20 ),OELTP, DELT 
1 1 BUT » 

1 NBV, 

1 ITRACE, 


135),XACTBV 

SVEC { 462),LF,IC,IIC,C( 462 ) , Y ( 20, 21 ) 
YV£C( 420 ) , YNL ( 80 ) , LF ( 20 ) , S { 2 1 , 22 ) , 
SPTRB, PNBRM, TSTART, TSTOP, 

SMALL, ' ZERO, DETT, TRACE, 

CONV, DPTRB, DTEMP, DN1, 

,TBV( 21),BV{ 20), 

, AUX { 10 ) ,CQBV( 20 ), 

M, TNORM,SN0RM,ANORM, TB VK , DOT, 

LIN, IQBVS, LITLR, nORMRD, 

NEQ, NCN, NWRITE, SKIP, 

LIMIT, IVARY, I A F ( 2 ) , 



£36 

837 

238 

239 
349 
£41 
34g 

243 

244 

245 

246 
£47 
24B 
249 
§50 
251 
25? 

253 

254 
?55 

256 

257 

258 

259 

260 
261 
262 
263 
§64 
§65 
266 
g67 
268 
269 
-SZS. 

271 

272 
?73 

274 

275 

276 

277 

278 

279 

280 
281 
28g 
?83 

284 

285 

286 

287 

288 

289 

290 

291 
392 
?93 
294 


1NVECT * IXTRA* 

1 JXTRA*KXTRA*N8UT* NL* N5PACE# NR CD* 

2NRSD*NRYD*I*II* IC0UNT * I N1 * 1 N3 * I NIKE Y* IPUNT* IQ* I S* I SUB* I TEMP* I TER* 

3 J* JJ* JQ* JSUB*K*KK*K0UNT* K0UT*KQ* KTEM * M I TER* NBVM* NC0L* NRANKC* 

4 NR0WS* QBV( 20 ) * XACTBV ( 21 ) *XACTI C( 2O>*IOBV!2O)*NR*IG0*NE*NI 
EQUIVALENCE ( YVEC* Y )* (S*SVEC ) * ( YVEC*C )* { I OUT* IB IG ) * ! DEL T* R0 IG ) 

34 IF (NBV 142* 42*919 

919 IF( NBV-1 >41*41*920 

920 IS=1 
NBVM«NBV-1 

35 DO 40 I=IS*NBVM 

IF ( TBV! 1+1 ] »T B V { I ) >37*40*40 

37 TEMP»TBV(1) 

TBV ( I )=TBV( 1+1 ) 

TBV < I+l ) *=T£MP 
TEMP=BV ( I ) 

B V ( I ) =B V { I + 1 ) 

B V ( 1 + 1 >»TEKP 
I TEMP=*QBV ( 1 > 

QBV ( I ) *QB V ( 1+1 ) 

QBV ( I + 1 ) = I TEMP 
ITEMP=XACTBV( I ) 

XACTBV ( I ) =X ACTBV ( 1 + 1 > 

XACTBV ( 1 + 1 i *1 TEMP 
I F ( 1-1)35*35*38 

38 I S= 1-1 
GO TO 35 

40 CONTINUE 

41 I F i TSTOP-TBVI NBV) >921*921*922 

921 TST0P=TBV ( NBV ) 

922 IF( TSTART»IBV( 1 > >42*42*923 

923 TSTART=TBV( 1 ) 

42 CONTINUE 

RETURN 

EM „ 

SUBROUTINE NTLIZE 
C 

INTEGER Q BV * XACT IC* SKIP* I BIG( 135)* XACTBV 

REAL RBIG(204>* LOWER*SVEC( 46 2 ) * LF* I C * I I C* C ( 462 ) * Y ( 20* 21 ) 

COMMON I IC ( 20 ) * IC< 20 ) * YVEC ( 420)*YNL( 80 > * LF ( 20 ) * S ( 2 1* 22 ) * 

1 ' DELT* SPTRB* PNORM* TST ART* TST0P* 

1 !)ET* SMALL* ZERO* DETT* TRACE* 

2 DZER0* C0NV* DPTRB* OTEMP* DN1 * 

3H*RN1*RN2*T* 

4TAVG*TL*TN* TBMT* TEM* TEMP* TTEM* TBV < 21)*BV( 20>* 

5 UPPER! 20 )* LOWER* 20 )* PTRBl 20) * AUX( 10 ) *CQBV(20 >* 

5 QIC! 20 )* YPTRBf 20 ) * DELTP* DELTM* TNORM* SNORM* ANORM* TBVK*DOT* 

1 I OUT* LIN* IQBVS* LITER* N0RMKD* 

1 NBV* NEQ* NCN* NWRIfE* SKIP* 

1 ITRACE* LIMIT* I VARY* IAF<2>* 

1NVECT* IXTRA* 

lJXTRA*KXTRA*NflUT*NL*NSPACE*NRCD* 

2NRSD* NR YO* I * 1 1* ICOUNT* INI * IN2* INI KEY* IPUNT* IQ* IS* I SUB* I TEMP* I TER* 

3 J* JJ* JQ* JSUB*K*KK* K0UNT* KOUT*KQ* KTFM*MITER*NBVM* NCOL/NRANKC* 

4 NROWS* QBV! 20 > * XACTBV! 21 ) *XACTICf 20 ) * I QBV ( 20 ) *NR* I GO* Nt*NI 
EQUIVALENCE ( YVEC*Y>* <S*SVEC>* (YVEC*C1* ( rflUT*IBlG>* < DEL T* RB I G > 
DELTP»1*1 *DELT 

DELTM = 0 *9 *DELT 

C#***INITIALIZE I C* D I C* AND XACTIC**** 



IY 


£95 


NVECT=NEQ+NCN 

296 

46 

DO 50 I-l/NVLCT 

297 


DIC(I ) = o. 

298 


ICC I )=IIC(I) 

£99 


IF(NBV)47#47#49 

300 

47 

I F ( XACTICt I ) ) 49# 48# 48 

301 

48 

X AC T I C ( I )=-l 

302 

49 

I F t XACTIC t I > + l 1924# 50# 50 

303 

924 

XACTICd) *= XAGT IC ( I ) +4 

304 

50 

CONTINUE 

305 


T NORM = 1. 

906 


SN0RM = PNORM 

30 7 


NR0WS=NBV+1 

3Q8 


IF ( NBV ) 61#61# 52 

309 

52 

DO 60 I=1#NVECT 

310 


•IF ( ABS< TBVI I ) -TSTART ) )55# 55# 61 

311 

55 

I-0 = QB V ( I ) 

312 


I F ( I0)57#57#56 

313 ' 

56 

I C ( IQ )=BV{ I ) 

314 

57 

I F { XACTBV I I ) ) 58# 58# 60 

316 

58 

I F ( IQ ) 60# 60# 59 

316 

59 

NROWS = NROWS -1 

317 


XACTBV ( I )=-l 

318 


X ACT I C ( I Q ) » X ACT I C ( IQ) - 4 

319 

60 

CONTINUE 

320 

61 

KOUNT = 0 

321 


DO 70 I=1#NVECT 

322 


I F ( XACTICf I ) ) 62# 63#*63 

323 

62 

KOUNT-KOUNT+1 

324 


GO TO 70 

325 

63 

IF ( A8S( PTRB ( I ) ) -DZERO 164# 70# 70 

326 

64 

PTRB [ I ) =SPTRB 

327 

70 

CONTINUE 

328 


ITER=0 

329 


M I TER*L I TER 

330 


IFCKOUNT-NVECT ) 926# 925# 925 

331 

925 

MITER=1 

332 

926 

NRANKC=NVECT-K0UNT+1 

333 


NC0L=NRANKC+1 

934 


RETURN 

_ 23 £_ 


■ENJD 


336 

337 
338' 
339 
349 

341 

342 

343 

344 
343 

346 

347 

348 

349 
•359 

351 

352 
35? 


SUBROUTINE ECHO IN 
C 

C ECHO I N 

INTEGER GBV#XACTIC# SKIP# I BIG ( 135) # XACTBV 

REAL RBIGI2Q4 !# LOWER . SVEC ( 462 ) # L F # I C # 1 1 C# C ( 462 ) # Y ( 20# 21 )' 

COMMON IIC, I20)#IC(H0)#YVEC{ 420>#YNL( 80 ) # LF < 20 ) # S ( 2 1 # 22 ) # 

1 DELT# SPTRB# PNORM# TSTART# T STOP# 

1 DEI# SMALL# ZERO# DETT# TRACE# 

2 DZERO# C0NV# DPTRB# DTEMP# DN1 # 

3H# RN1#RN2# T# 

4‘TAVG#TL#TN#TBMT#TEM#TEMP J ,TTEM#TBV< 21 ) # B V < 20 1# 

5 UPPER(20)#LOWER(20)#PTRB(?0)#AUX( 1 0 ! # COB V ( 20 ) # 

5 OIC(2O)#YPTRB(2O)#DELTPjDELTM#TNORM#SNflRM^AN0RM,TBVK#D0T> 

1 I0UT#' LIN# IQBVS# LITER# NORMHD# 

1 NBV# NFQ# NCN# NWRITE# SKIP# 

1 ITRACE# LIMIT# IVARY# IAF(2)# 

1NVECT# IXTRA# 

1JXTRA#KXTRA#^8UT#NL#NSPACE#NRCD# 
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35* 

355 

356 

357 

358 

359 

360 

361 

362 

363 
36* 

365 

366 

367 

368 

369 

370 

371 

372 
379 
37* 
975 

376 

377 

378 

379 

380 

381 

382 

383 
38* 

385 

386 

387 

388 

389 

390 

391 

392 
399 
39* 

395 

396 

397 

398 

399 
*00 
*01 
*02 
*03 
*0* 
*05 
*06 
*07 
*08 
*09 

*11 

*12 


2NRS D* NR YD/ I * I 1/ I COUNT * INI * I N2 * I NT KEY* IPUNT* IQ* IS* ISUB* I TEMP* I TER, 

3 J* JJ* JU* JSUB* K*KK*K0UNT* K0UT * K0* KTEM* M ITER* NBVM* NCOL*NRANKC* 

* NR0WS*QBV( 20 } * XACT BV <21 ) * X AC T I C ( 20 ) * I8BV(20t*NR* IG0*NE*NI 
EQUIVALENCE ( YVEC*Y)* (S*SVECI* ( YVEC*C>* { I OUT* IBIG)* ( DEL T * RB I G > 

C 

C*#**ECHS THE INPUT DATA*»x* 

WRITE (6*350) 

I E ( NB V ) 73*73* 927 

927 DO 72 1 = 1* NBV 

WRITE <6,»*00) I*TBV ( I >*BV< I )*QBV( I )*XACTBV< I ) 

IF (QBVl l ) >928*72*72 

928 J» -QBV ( I 1/100 

K = -QBVl I ) - j-noo - 1 + J 

WRITE (6***0) ( CQBV ( M ) * IQBV (H ) * M* M^J* K ) 

72 CONTINUE 

73 WRITE (6*450) 

DO 7* I =1 * NVE CT 

7* WRITE (6*500) I * I C ( I ) * X ACT I C i I ) t PTRB ( I ) * UPPER ( I ) * L0 WER < I ) 

WRITE (6*307 ) I0UT* LIN* IQ BVS* LITER* N0RMRD * NBV * NEQ* NCN* N WR I TE* SK 1 P* 
1DELT*SPTRB* PN0RM* TSTART* TST0P* NRANKC 
360 FORMAT <*H AUX* /* ( IX* 8E 15 • 7 ) ) 

WRITE (6*360) < AUX ( I ) * I =1 * 1 0 ) 

RETURN 

307 FORMAT (//1OX*4HI0UT* 16* 10X*3HLIN* 17* 5HIQSVS* 15* 10X* 5HLITER* lb* 10X 
1* 6HN0RMRD* I */ / 1 2X * 3HNBV * 17* 10X* 3HNEQ* 17* 10X* 3HNCN* 17* 10X* 6HNW *1 TE* 
21** 10X* *HSKIP* I6//1X* *HDELT*E16t6*5X* 5HSPTRB* E 1 5 . 6* 5X* 5HPN0RM* £ 1 5 ■ 
36* 5X* 6H TSTART* El* .6*5X* 5HTST0P* El 5 * 6/// 1 OH GENERATED/10X* 6H NRANKC* 
*1*) 

350 FORMA T( IX// 16X*3HTBV* 17X* 3HB-V* I3X* 3HQBV* *X*5HXCTBV> 

400 F0RMAT ( 1 X* 15* 2 I 5X* El 5 . 7 ) * 2 18 ) 

4*0 FORMAT <30X*E15»7*3H*Y ( * 14* 1H )*20X* 15 ) 

4 50 F0RMAT ( 1X//3X*3HY( ) * 8 X* 2H I C * 1 3X* 5HXCT IC* 1 OX* 4 HP TR B* 1 5 X* 

1 5HUPPER* 15X* 5HL0WER ) 

500 FORMAT! *X* 12* 5X> E 15 * 7*5X* 13* 4 ( 5X* E15 . 7 ) ) 

- 

SUBROUTINE ESTABL 
C 

INTEGER QBV * X ACT I C* SKIP* IBIG(135)*XACTRV 

REAL RBIG120* )* LOWER* SVEC ( *6 g ) * LF* IC* I IC* C ( *62 ) * Y ( 2 0* 21 ) 

COMMON IIC(20)*IC(20)*YVECl *20)*YNL( 80 1 * LF ( 20 ) * S ( 2 1* 22 ) * 

1 DELT* SPTRB* PN0RM* TSTART* TSTOP* 

1 DET* SMALL* . ZERO* DETT* TRACE* 

2 DZERO* C0NV* DPTRB* DTEMP* DN1* 

3H-» RN 1 * RN2* T * 

4TAVG*TL*TN*TBMT*TEM*TEMP*TTEM*TBV( 21)*BV( 20)* 

5 UPPER ( 20 ).* LOWER ( 20 )* PTRB ( 20 ) * AUX ( 10 ) *CQBV( 20 )* 

5 DI C ( 20 ) * YPTRB ( 20 ) * DEL TP* DELT M* TN0RM* SNORM* A NORM* TBVK*DOT * 

1 I OUT* LIN* IQBVS* LITER* NORMKD* 

1 NBV* NEQ* NCN* NWRITE* SKIP* 

1 ITRACE* LIMIT* IV ARY* IAF(2)* 

1NVECT* IXTRA* 

1 JXTRA*KXTRA*N0UT* NL'*N SPACE* NRCD* 

2NRSD* NR YO* 1 * I I* I COUNT * I N 1 * I N2* I NTKE Y* IPUNT* IQ* IS* I SUB* I TEMP* ITER* 

3 J* JJ* JQ* JSUB* K* KK* KOUN T*KOUT*KQ*KT EM* MI TER* NBVM*NCOL* NRANKC* 

* NROWS* QBV ( 20 ) * XACTBV ( 21 ) * X ACT I C { 20 ) * 1 68V ( 20 ) * NR* I GO* NE* N I 
EQUIVALENCE (YVEC*Y)* (S*SVECI* (YVEC*C)* < I 3U T* I B lG ) * ( D EL T* RB IG ) 

C****ESTABL ISH INITIAL CONDITIONS IN YVEC*#*V 
79 I TER=IT£R+1 
TRACE * 0. 
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413 
if 3.4 

415 

416 
517 
5-1 B 
519 

420 

421 

422 

423 

424 

425 

426 

427 

428 
429 
539 
431 
532 

-&a& 

434 

435 
536 
4 3.7 
538 

439 

440 
451 

442 

443 

444 

445 
-546 
447 
4.4 8 
459 

450 

451 

452 

453 

454 

455 

456 

457 

458 
.459 
46.0 
461 
562 
463 
465 
465 
566 

467 

468 

469 

470 
571 


JSUB=-NRYD 
JQ = 0 

DO 100 J* 1 / NR ANKC 
JSUB=JSUB+NRYD 
DO 80 I *1 / N VE'CT 
IZ = JSUB 4 1 
80 YVE C ( IZ)=IC.< I ) 

IF ( JSUB ) 82/ 100/ 82 
82 IS=JQ+J-1 

DO 90 I = I S/ NVECT 
IFI XACTICU ) ) 90/ 92/ 92 
90 JO = JQ 4 1 

92 JJ = J - 1 4 JQ 

I Z* JSUB 4 J J 

Y VEC ( IZ )=YVEC( IZ)*PTRB( J J ) 

IF t ABSI YVEC ( IZ )-DPTRB ) >.93/9 4/ 94 

93 YVEC I IZ > = 1. 

94 YPTRBl JJ)*YVEC( 12) 

100 CONTINUE 

RETURN 



SUBROUTINE NGRATE 
C 

INTEGER QBV/XACTIC/SKIP/ I BIG ( 1351/XACTRV 

REAL RBIG ( 204 I / LOWER/ SVEC ( 462 ) / L F / IC/ I I C/ C ( 462 } , Y ( 20/ 2 1 ) 

COMMON IIC(20)/1CI20)/YVEC( 42Q)/YNL< 80 > / LF ( 20 J / S ( 2 1 / 22 ) / 

1 DELT/ SPTRB/ PN0RM/ TSTARl / TSTOP/ 

1 DEI/. SMALL/ ZERO/ DETT/ TRACE/ 

2 DZERO/ - CON V/ DPTRB/ DTEMP/ DN 1 / 

•3H/RN1/RN2/T/ 

. 4TAVG/TL/TN/TBMT/TEM/TEMP/TTEM/TBV( 21)/ B V I 20)/ 

5 UPPER! 20 ) / LOWER! 20 > / PRT6 ( 20 ) / AUX { 10 ) /CQBVI 20 ) / 

5 DIC(2O)/YPTRB(2O)/DELTP/DELTM/TN0RM/SN0RM/ANORM/TBVK/D0T/ 

1 I OUT/ LIN/ IQBVS/ LITER/ N0RMKD / 

1 NBV/ NEQ/ NCN/ NWRIT-E/ SKIP/ 

1 ITRACE/ LIMIT/ I VARY/ IAF(2)/ 

INVECT/iXTHA/ 

1 JXTRA/ KXTRA/ N0UT/ NL/NSPACE/ NRCD/ 

2NRSD/ NR-YD> I / 1 1/ I COUNT / I N 1 z I N2/ I NTKE Y/ IPUNT/ 10/ I S/ I SUB/ I TEMP / I TER/ 

3 J/ JJ/ JQ/ JSUB/K/KK/ KOUNT/ K0UT/KQ/ KTEM/MITER/NBVM/NCOL/NRANKC/ 

4 NROWS/QBV( 20 )/ XACTBV ( 21 ) /XACTIC( 20 )/ IQBV { 20) /NR/ IGO/NE/NI 
EQUIVALENCE ( YVEC/Y )/ (S/SVEC) / ( YVEC/C )/ ( I0UT/ 1 8 I G ) / ( DEL T/ RB I G > 

c****inte.gration**** 

102 K0UT“0 

I'F< I.0UT-2 ) 104/ 104/ 103 ' 

104 IF( ITER-MITER >929/103/ 103 

103 K0UT=1 

929 IFI 10UT-4 >931/930/930 

930 WRITEI6/630) 

931 NI=NR0WS41 
N E *‘l 

JSUB*-NRSD 

DO 105 1*1/ NCOL 

JSUB=JSUB4NRSD 

105 SVECI JSUB41 }«l . 

TL=TSTART 
ICOONT=NWRI TC ‘ 

IF( ABSITBVI 1 ) -TST ART) -ZERO) 106/ J 06/ 107 

106 IC0UNT»0 
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£72 

473 

474 
£75 
£76 
£77 

478 

479 

480 

481 
£82 

483 

484 

485 

486 
£87 
£88 
489 
£99 
£91 
£92 
£93 
494 
£95 

496 

497 
£98 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 
611 

512 

513 

514 

515 

516 

517 

518 

519 

520 

521 
622 

523 

524 

525 

526 

527 

528 

529 

530 


107 IF< ITER“MlTER)936/935/935 

935 WR I TE <6/590 I 

936 WRITE (6/600 ) I TER/MI TER 
I NTKE Y=0 

INI = NRANKC 

932 IF( MITER-lTER)933/933/934 

933 I N1 = 1 

934 K S 1 

TBVK * TBV ( K > 

110 TEM«0» 

RETURN 

600 FORMAT f 20X/4HITER/ I4/4H OF/ 14/7/ 10X/4HTIME/ 14X/8H SOLUTION) 

630 F 0RMA T ( 1 X / l'OX/ 3 9H RE PE A TED TIMES ARE PARTICULAR SOLUTIONS//) 

590 FOR MA T ( 1 H 1 ) 

EJiC 

SUBROUTINE STORES 
C 

c STORES 

INTEGER QBV/XACTIC/SKIP/ IBIGi 1351/XACTBV 

REAL RBIGI204 )jLOWER/SVEC( 46 2 ) j LF/ I C / I IC/ C ( 462 ) / Y < 20/ 21 > 


COMMON 

I IC ( 20 ) / I'CI 201/YVEC ( 420) 

/ YNL ( 

80)/LF(20)/S(21/22) 

1 

DELT / 

SPTRB/ 

PN0RM/ 

TST ART/ TST0P/ 

1 

DET/ 

SMALL/ 

ZERO/ 

DETT/ TRAGE/ 

2 

3H/RN1/RN2/T/ 

PZERO/ 

CONV/ 

DPTRB/ 

DTEMP/ DN 1 / 


4TAVG/ TL/TN/ TBMT/TEM/TEMP/ TTEM/TBV ( 21')/BV( 20)/ 

5 UPPER ( 20 ) / LOWER ( 20 ) / PTRB ( 2 0 ) / AUX ( 10 ) /CQBV( 20 ) / 

5 DIC( 20 )/ YPTRB ( 20 ) /DELTP/ OELTM/ TNORM/SNORMj ANORM/TBVK/DOT/ 

1 I OUT/ LIN/ IQRVS/ LITER/ NORMRD/ 

1 NBV/ NEO/ NCN/ NWRITE/ SKIP/ 

1 ITRACE/ LIMIT/ IVARY/ IAF(2)/ 

1NVECT/IXTRA/ 

1 JXTRA/KXTRA/NOUT/ NL/NSPACE/ NRCD/ 

2NRSD/NRYD/ 1/ II/ 1C0UNT/ INI/ IN2/ INTKEY/ IPUNT/ IQ/ I S/ I SUB/ I TEMP/ I TER/ 

3 J/ JJ/ JQ/ JSUB/K/KK/ KSUNT/ KOUT/KQ/ KTEM/MITER/NBVM/NCOL/NRANKC/ 

4 NR0WS/QBV{2OI/XACTBV(21)jXACTIC(2O)/IBflV(2O)/NR/IG0/NE/NI 
EQUIVALENCE ( YVEC/ Y ) / (S/SVEC) / ( YVEC/C )/ ! I OUT/ IBIG )/ (DELT/ RBIG ) 

134 IF (XACTBV(K)) 152/135/135 

135 10 * QBV(K) 

IF ( XACTBV (K ) ) 139/ 138/ 139 


139 

Nl=NI-l 

nr=ni 

GO TO 140 


138 

NE = NE + 1 
NR = NE 


140 

I SUB=-NRYD 
JSUB=-NRSD 
IF { IQ) 143/ 143/ 144 


143 

J- -IQ/100 



IS « J - (IQ + J*IOO ) - 

144 

DO 150 I=1/NRANKC 


1 SUB » I SUB + 

NR YD 


JSUB = JSUB + 
I F ( IQ )14b/ 145/ 147 

NRSD 

147 

IZ * IQ + ISuB 
DOT = Y VEC ( IZ1 
GO TO 148 


145 

DOT = 0. 



IF <IS-J)149/950/9bO 



§31 

950 

00 146 M*=J#IS 

532 


JJ ■ IQBV(M) 

533 


IZ * JJ + 1SU8 

534 

146 

DOT = DOT + YVEO 1Z)*CQBV('M! 

535 

148 

IZ = NR + JSUB 

536 


SVEC(IZ) » DOT 

§37 

149 

IF ( I-NRANKC ) 150# 951# 150 

533 

951 

IZ = NR+NRSD+JSUB 

539 • 


SVEC< IZ ) * BV (K ) 

§40 

150 

CONTINUE 

541 

151 

I F ( IQ ) 952/952# 152 

542 

952 

I F ( IS “J)9 53 #152/152 

543 

953 

CONTINUE 

544 


CALL NLBC 

545 

152 

K *» K +• t 

546 


IFf K-NBV ) 154# 154# 153 

547 

153 

IF) ITER-MITER ) 170# 154# 154 

548 

154 

IF) K-NBV ) 954# 954# 955 

•54 9 

954 

TBVK^TBV ( K ) 

550 


60 TO 956 

551 

955 

TBVK=TST0P+2»*DELT 

552 

956 

I F < AB S { TB VK "*T B V ( K - 1 ) ) ) 134/134/957 

553 

957 

TEM=TEM+TBMT 

554 


TEM=TEM+T8MT 

555 


ttem=tl+delt-tem 

556 


IF ( TEM-DELTM) 155/958# 958 

557 

958 

TTEM=TL+2»*DELT-TEM 

558 


• KTEM=KTEM + 1 

559 


I F ( KTEM-1 >960/959/960 

560 

959 

ICQUNT=ICBUNT + 1 • 

561 

960 

I F < I COUNT -N WRITE ) 155# 961 # 15 5 

562 

961 

I COUNT = 0 

563 

155 

INTKEY=0 

564 


X F { ABS( TTEM-TBVK) •) 110/110/157 

565 

157 

IF! TTEM-TBVK >158#158#112 

566 

158 

h=ttem-tl 

567 


I=-l 

568 

110 

RETURN 

569 

170 

I = 0 

570 


RETURN 

571 

112 

J = 0 

572 


RETURN 

Jaza 


-ZfciiJ . 

574 


SUBROUTINE NLBC 

575 

C 


576 


REAL LF/DY) 20 )#U( 20 ) 

577 


COMMON F1C( 40)# P( 20/21 ) #0 80) /LF( 20)# 3(21/22)# 

578 


1 R01 ( 9 ) # TRACE#RD2(5 )#H# R03( 2 ) # T^RD4 (7 ># TBV( 21 )#BV( 80 ) # 

579 


2 AUX ( 10 ) # CQ ( 67 ) # 

580 


3 I 01 { 16) # I PART# JSUB# 1 02 ( 2 ) # NL# I D3 ( ?9 ) # MR ANKC# I 04 ( 6 2 ) # I Q ( 24 ) 

581 


EQUIVALENCE (TL#RD4(2))#(J#ID3(17))#(K#ID3(21))#(NR#IQ(21)) 

582 


EQUIVALENCE ( LF# DY ) # ( P# U ) 

583 


DO 80' I CL » 1 /NRANXC 

584 


ASM = 0 • 

585 


DO 79 IDX * 1/4 

586 


IZ=5-IDX 

587 

79 

ASM-ASM-U t IZ > *P ( IDX# I CL 1 

588 

80 

SINR/ICL) » ASM 

589 


S ( NR# NRANKC + 1 )*8V1K 1- S IN ( U ( 5 ) *TL > -U <1 )-*U ( 4 ) -U ( 2 ) *U ( 3 ) 



590 


RETURN 

591 


END 

592 


SUBR0UT INE NEWIC 

593 

C 


594 

C 

NEWIC 

595 


INTEGER QBV/XACTIC/SKIP/IBIS< 135) /XACTBV 


596 

REAL RBIG(204)/L0WER/SVEC( 462)/LF/IC 

/IIC/C( 462 ) / Y ( 20/ 2 

597 

COMMON IIC(20)/IC(20)/ YVECf 420)/YNL< 

80 ) 

/LF ( 20 ) / S ( 2 1/ £2 ) 

59S 

1 CELT / SPTRB/ PN0RM/ 


TSTA'RT/ TSTOP/ 

599 

1 DET/ SMALL/ ZERO/ 


DETT/ TRACE/ 

£0 0 

2 DZtRO/ C0NV/ DPTRB/ 


DTEMP/ DN 1 / 

601 

3H / RN 1 / RN2/ T / 



602 

4TAVG/TL/TN/TBMT/TEM/TEMP/TTEM/TBV( 21 ) / BV { 20 

)/ 


603 

5 UPPER ( 2 0 ), LOWER ( 20 ) / PTRB ( 20 > / AUX ( 1 0 ) / CQBV ( 20 

)/ 


604 

5 DICI 20 ) j YPTRB( 20 1 / DEL TP/ DEL T M/ TN8RM/ SN BRMj AN 

ORM 

/TBVK/DOT/ 

§05 

1 I OUT/ LIN/ IQBVSj 


LITER/ NO RMRD / 

606 

1 NBV/ NEQ/ NCN/ 


NWRITE/ SKIP/ 

607 

1 1 TRACE/ LIMIT/ IVARY/ 

I AF ( 2 ) / 

§08 

1NVECT/ IXTRA/ 



§09 

1JXTRA/KXTRA/N8UT/NL/NSPACE/NRCD, 



610 

2NRSD/ NR YD/ I / I 1/ I C BUNT/ IN1/INH/INTKEY/IPUNT/IQ 

/IS 

/ ISUBj ITEMP/ iter 

611 

3 J/JJ/JOj JSUB/K/KK/ (COUNT/ KOUT/KQ, KTEM / MITER / NBVM/NCSL/NRANKC/ 

612 

4 NR0WS/ QB V < 20 )/ XACTBV (21)/XACTIC(20)/I0BV(20) 

/NR 

/ IG0/ NE/NI 

613 

EQUIVALENCE ( YVEC/Y )/ (S/SVEC) / ( YVEC/ C ) / ( 1 0U T/ 

IBIG)/ (DELT/RBlG) 

614 

C*#**CALCULATE NEW INITIAL COND I TIONS**** 



615 

IF "( J) 970/ 188/ 188 



616 

188 KQ=0 



617 

AN0RM = 0» 



61B 

DOT » 0* 



619 

JSUB=NRANKC*NRSD 



620 

DO 210 J=1/NVECT 



621 

IF(XACTIC(J))19 0/195/ 195 



622 

190 KQ=KQ+1 



623 

G0 TO 210 



624 

195 JJ*J-KQ+1 



625 

C****CALCULATE NEW J-TH ELEMENT OF INITIAL CONDITION 

VECTOR**** 

626 

IZ - JJ + JSUB 



627 

DTEMP* SVEC ! IZ) 



628 

DTEMP = YPTRB ( J )*DTEMP + IC(JJ*(1» -DTEM.P) 



629 

TDIC =. DTEMP - IC ( J ) 



630 

I F ( A8S( PN0RM) -ZERO) 205/ 198/198 



631 

198 DTEMP=DIC1 J ) *0 • 5 



632 

IF(TDIC*DIC(J))199/ 205/205 



633 

199 IF ( ABS ( TDIC ) -ABS 1 DTEMP M 205/205/200 



634 

200 TDIC*-DTEMP 



635 

205 DOT * DOT + TDIC*DIC(J) 



636 

ANBRM = ANORM + TDIC**2 



637 

DIC(J) = Tv'DIC 



638 

210 CONTINUE 



639 

ANORM = SQRHANORM) 



640 

DOT * DOT/( ANORM+TNORM) 



641 

212 IF( PNORH) 970/ 970/213 



642 

213 CONTINUE 



643 

J «= «*1 



644 

RETURN 



645 

9 70 I F < ANORM-SNQRM) 97 1/963/ 96 3 



646 

971 I F { N0RMRD J963/963/962 



647 

962 SN8RM=AN0RM 



648 

963 'TN0RM = AN0RM 





649 

£51? 

651 

652 
659 

655 

656 

65 7 

658 

659 

660 
661 

66 2 
■ 663 

664 

665 

666 

668 

669 

67Q 

671 

672 

673 

674 

675 

676 

677 

678 

679 

680 
681 
682 

683 

684 

685 

686 

687 

688 

689 

690 

691 

692 

693 

694 

695 

696 

697 

698 

699 

700 

701 

702 

703 

704 

705 

JLM 

707 


00 215 J = 1 * N VECT 
I F ( XAC7 I C ( J ) > 215/ 96 4/ 96 4 

964 IC( J)»IC( J)+DIC(J) 

IF (XACTIC1 J)-l 1214/965/ 965 

965 IF ( IC ( J > -UPPER ( J ) >967*967*966 

966 I C ( J ) “UPPER ( J ) 

9 67 IF(IC(J}-L0WER(J) > 968 j 214/214 

968 1 C ( J ) =L0WER ( J ) 

214 I F t SKIP >215*969*215 

969 HCtJ)-IC(J) 

215 CONTINUE 

DO 220 1=2* NRANKC 

I Z= I + JSUB 

IF( ABS(SVECdZ) )«C0NV >220*216*216 

216 X F ( ITER-LITER >79* 220* 220 
220 CONTINUE 

MITER=ITER+1 
79 RETURN 

END 

SUBROUTINE MBDIFY 
C 


INTEGER QBV*XACTBV*SKIP* I B I G < 135) 

REAL RBIG (204 )* LOWER* SVEC ( 462 ) * LF* IC* I IC* C ( 462 ) * Y ( 20* 2 1 ) 

COMMON I IC ( 20) * IC ( 20 )* YVEC( 420>*YNL< 80 ) * LF ( 2 0 ) * S ( 2 1 * 22 ) * 


1 

DELT* 

SPTRB* 

PNORM* TSTART* 

TSTOP* 

1 

DET* 

SMALL* 

ZERO* DETT* 

TRACE* 

2 

DZERO* 

CONV* 

DPTRB* DTEMP* 

DN1* 

3H*RN1*RN2*T* 




4TAVG*TL*TN*TBMT*TEM*TEMP* 

TTEM* TBV ( 2 1 ) * BV ( 20)* 


5 

UPPER! 20 >* LOWER ( 20 ) * PTRB { 20 ) * AUX ( 10 ) 

* CQB V ( 20 ) * 


5 

DIC( 20 )* YPTREM 20 ) *DELTP* 

DELT M* TNORM* 

SN0RM* AN0RM* TB VK* DOT * 

1 

TOUT* 

LIN* 

IQBVS* LITER* 

N6RMK0* 

1 

NBV* 

NEQ* 

NCN* NWRITE* 

SKIP* 

1 

ITRACE* 

LIMIT* 

IVARY* I A F ( 2 ) * 


1NVECT* IXTRA* 

. 1JXTRA*KXTRA*N0UT*NL*NSPACE*NRCP# 

2NRSD*NRYD* I* I I* I C OUNT * I Nl * I N2* I NT KE Y* IPUNT* 10* IS* I SUB * I TEMP* I TE R* 

3 J*JJ*jQ*JSUB*K*KK* KOUNT* K0UT*KQ* KTEM* MITER* NBVM* NCOL* NRANKC/ 

4 NROWS* QBV ( 20 ) * XACTBV ( 20 ) * K I C ( 2 1 ) * I QBV ( 20 ) * NR* I GO* NE * N I 

EQUIVALENCE { YVEC* Y ) * { S* SVEC ) * ( YVEC* C ) * ! I OUT* IB IG ) * < DELT* RBIG ) 

IF (AN0RM-PN0RM )2*4*4 

2 IF ( DOT >4*4*3 

3 RETURN 

4 C0EF=1 • 

IF < AN0RM-PNORM ) 6*6*5 

5 COEF = PN0RM/AN0RM 

6 ITERM = -PNORM/AN0RM - 1. 

I F ( ITERM+10 >7*8*8 

7 I TERM*- 10 

8 IF{D0T)9* 10* 10 

9 COEF = COEFM 1 • +DOT ¥ 1 0 . I TERM ) 

10 DO 20 I=1*NVECT 

I F < K I C ( I ) >20*20* 11 

11 DIC(I ) = CBEF*DIC( I > 

20 CONTINUE 

ANORM = ANORM^COEF 
RETURN 

END 


SUBROUTINE GJRWLS 





708 

C 


709 


INTEGER I(21)*J(22)*J2C(20)*IJTV(20> 

710 


COMMON ADI (40 >*B( 21*22) *AD2(5 8)* A (21* 22)* AD3( 8 ) * SE T * T RA CE* CHE CK* AD 

711 

14(1 93 )*I0UT *101(21 ) *NRB*NRA* I 02 ( 7 ) * I PUNT* I D3 ( 1 8 ) * N* NR0WA*.ID4 ( 60 ) * 

712 


2 NEMAX* 105(24) 

713 


EQUIVALENCE < AD2< 39 )* IJTV > 

714 


NTOTAL N f 1 ■ 

715 


lst-op»nemax 

716 


NR0W=NR0WA 

717 


nend=nrowa' 

718 


IF (NROtiA-N TOTAL) 20*30*30 

719 

20 

NEND=NT8TAL ' 

7S0 

30 

00 100 L= 1 * NE ND 

731 


IF (.L-NT0T AL ) 4 0* 40* 100 

722 

40 

J ( L ) a L 

733 

100 

I ( L ) =L 

724 


00 650 K=1*N 

725 


C = 0«0 

726 


D 0 200 L = K* LSTOP 

727 


IL = ' I ! L) 

728 


D S 200 M = K* N 

729 


JM=J:(M) 

730 


ABS A= ABS { A ( IL* JM) ) 

731 


I F { ABSA-C ) 200* 200* 190 

732 

190 

KL = L 

733 


KM ® M 

734 


C = ABSA 

735 

200 

CONTINUE 

736 


KI » I ( KL ) 

737 


K J = J ( KM ) 

738 


I (-KL ) = I ( K ) 

739 


I < K ) = KI 

740 


J(Ktl) = J(K ) 

741 


J ( K ) = KJ 

742 


AKK = A ( K I * K J ) 

743 


A(KI*.KJ)-1» 

744 


DET * DET* AKK 

745 


IF( NEMAX-N ) 26 0* 260* 250 

746 

250 

I F ( K-N ) 260* 1400*1400 

747 

260 

IF ( C-CHECK ) 345* 34 5* 350 

748 

345 

CONTINUE 

749 


I PUNT a 2 ' 

750 


RETURN 

751 

350 

AD2 IK + 18 ) » A KK 

752 


IJTV(K)= 1Q0*KI+KJ 

753 

370 

AKK * 1 • EO/ AKK 

754 


KP1=K+1 

755 


DO 40Q M= KP 1 * NT 3T AL 

756 


UM=J(M) • 

757 

400 

A(KI*JM)=AKK*A(KI*JM) 

758 


DO 600 L=1*NR0W 

759 


IF ( L"K ) 420* 600* 420 

760 

420 

IL*I( L) 

761 


AKK = A ( I L* K J ) 

762 


A ( IL*KJ 1=0. 

763 


DO 500 M=KP1*NT0TAL 

764 


JM = J ( M ) 

765 

'500 

A ( I L* JM ) = A ( 1L* JM ) - AKK ♦A (K I #. V JM ) 

76 6 

60 0 

CONTINUE 



25 


767 

768 

769 

770 

771 

772 

773 

774 

775 
7 76 

777 

778 

779 
789 

781 

782 

783 

784 

785 

786 

787 

788 

789 

790 

791 
7?e 

793 

794 

795 

796 
79? 

798 

799 
SQ0 
801 
802 

809 

804 

805 

806 

807 

808 
8Q9 

810 
811 
812 
613 

814 

815 

816 

817 

818 

819 

820 
821 
822 

823 

824 

825 


IF ( K-NEMAX I 650/ 610/ 650 
610 IF(N«K)650/650/620 

620 00 621 L = 1/N 
JZC(L) *0 

DO 621 M= 1 / NTOT AL 

621 B(L/M) - 0. 

DO 627 L=KP 1 / N 
IRAT * J(LJ 
JZCIIRAT) = 1 
00 626 M*KP1/NT0TAL 
JCA = J(M) 

IF (IRAT-JCA) 622/ 622/ 626 

622 TEMP = Ot 

D 0 623 KR = KP1/NR6WA 

623 TEMP = TEMP + A (KR/ IRAT )*A(KR/JCA > 

BtIRAT/JCA) = TEMP 

IF (IRAT-JCA.) 624/626/626 

624 IF (JCA— N) 625/625/626 

625 B ( JCA / IRAT) « TEMP 

626 C0NTINUE 

627 CONTINUE 
IRA * K 

00 630 L * 1/N 

IF ( JZC ( L ) ) 628/630/628 

628 IRA = IRA + 1 

00 623 M =*= KP1/NT0TAL 
JCA = J ( M ) 

629 AdRA/JCA) = B(L/JCA) 

630 CONTI NUE 
LST0P = N 
NR0W = N 

650 CONTINUE 

DO 1000 K = 1 / N 
K I = I ( K ) 

KJ * J(K) 

IF(KI-KJ)660/800/660 
660 DO 700 M-l/NTOTAL 
AKK=A ( KI/M ) 

AIKI/M-) = A (KJ/ H) 

700 A ( K J/ M ) = AKK 
800 DO 900 L ~ K/N 

1 F ( I ( L ) '* J ( K ) 190 0/ 95 0/900 
900 CONTINUE 

950 I I L ) * I ( K ) 

1000 I (K ) » L 
1400 CONTINUE 
1450 I P U N T * 1 
RETURN 

END 

SUBROUTINE ONETEN 

INTEGER QBV/XACTIC/ SKIP/ 1BIG < 135 ) /XACTBV 

REAL RBIG( 204 )/ LOWER/ SVEC( 462 ) / LF/ IC/ I IC/ C ( 462 ) / Y ( 20/ 21 ) 


COMMON I IC ( 20 ) / I C( 20 ) / YVEC ( 420)/YNL< 80 I / LF ( 20 ) / S < 2 1 / 22 > / 


1 

DtLT/ 

SPTRB/ 

PNORM/ 

TSTAKTz 

T STOP/ 

1 

DLT / 

SMALL/ 

ZERO/ 

OETT/ 

TRACE/ 

2 

OZLRO/ 

C 0 N V / 

DPTRB/ 

DTEMP/ 

DN1/ 


3H/RN1/RN2/T' 

4TAVG/TL/TN/T0MT/TEM/TEHP/ TTEM/TBV( 2l)/BV( 20)/ 



26 


826 

827 

828 

829 

830 

831 

832 
838 

834 

835 

836 

837 

838 

839 

840 

841 

842 

843 

844 

845 

846 

847 

848 

849 

850 

851 

852 

853 

854 

855 

856 

857 

858 

859 

860 
861 
862 

863 

864 

865 

866 

867 

868 

869 

870 

871 

872 

873 

874 

875 

876 

877 

878 

879 

880 
881 
882 
B83 
884 


5 UPPER ( 20 > /LOWER <20 1/ PR TB (201* AUX( 1 0 ) /CQBV( 20)/ 

5 DIC (20)/ YPTRB<20 > / DELTP/ DELTM/ TN0RM/ SN0RM/ AN0RM/ TBVK/D6T/ 

1 1 0UT/ UN/ IQBVS/ LITER/ N0RHRD/ 

1 NBV/ NEQ, NCN/ NWRITE/ SKIP/ 

1 ITRACE/ LIMIT/ I VARY/ IAF(2)/ 

1NVECT/ IXTRA/ 

1 JXTRA/KXTRA/N0UT/NL/NSPACE/ NR CD/ 

2NRSD/NRYD / 1/ 1 1/ IC9UNTUN1 /IN2 / INTKEV/ IPUNT/ IQ/ IS/ ISU8 / ITEMP / I TER/ 

3 J/ JJ/ JQ/ JSUB/K/KK/K0UNT/ K0UT / K Q/ KTEM/ M I T ER / N B V M/ NC0L / NRANK C/ 

4 NR0WS/QBV(2O)/XACTBV(21)/XACTIC(2O)/J@BV(2O)/NR/IG0/NE/NI 
EQUIVALENCE ( YVEC/ Y ) / ( S/SVEC ) / ( YVEC/C )/ ( I0UT/ IBIG )/ 1 DELT/ RBIG ) 

G0 TO ( 111/ 105/ 105/ 105 >/ IVARY 

105 IF (J) 106/112/106 

106 IF (I) 118/110/110 

110 TEN- = 0. 

KTEM=0 

IF (TL-TST0P ) 937/937/135 

937 IF < IC0UNT-NWRITE > 112/938/ 938 

938 I G0» 1 

109 IF(K8UT-1 1939/940/940 
9 39 IF ( ABS( TL-TST ART) ) 940/ 940/941 

941 GO T 8 ( 111/134)/ IG0 
940 N0UT*NEQ 

lF( ABS( TL-TSTART) >942/942/943 

942 N0UT»NVECT 

943 CONTINUE 
IVARY = 1 
RETURN 

111 I C0UN T = 0 

112 TBMT = TBVK - TL 

IF( TBMT»OELTP 1115/ 115/ 114 
115 INTKEY=1 
H«TBMT 
G0 T 0 118 
114 H=DELT 
118 T = TL + H 
CALL DOIT 

tl=tl+h 

I F ( lNTKEY»l >131/130/131 
131 IC0UNTUC0UNT + 1 
G0 T 0 110 
130 I GO * 2 
G0 TO 109 

134 IVARY = ? 

RETURN 

135 IVARY = 3 
RETURN 

EJVD 

SUBROUTINE D0IT 
C 

REAL YK(201/YP(20)/A(3)/TN(4)/B.(3)/LF 

COMMON DDK 40)/ Y( 420)/ C ( 80 ) / LF { 20 ) / SVEC < 462 ) / 

1 RD1 ( 15 ) /HH/RD2 ( 4 ) / T/RD3( 183 )/ 

2 10UT/ LIN/ IDK 4 ) /NEQ/ 102 ( 8 ) / MVECT/ I PART/ JSUB/ 

3 1 03 ( 2 ) / NL / I D4 ( 3 ) /.NPYD/ I D5 ( 3 ) / IN 1 / T 06 ( 1 08 ) 

0 AT A A/1./2./2*/ 

I P ART*0 

iSUBa-NRYO 

1 G0=O 



885 

1 

IPART*=IPART + l 

886 


T I = T 

887 


IF< IPART-IN1 >3*3*2 

888 

2 

RETURN 

889 

3 

ISUB=ISUB+NRY0 

890 


H=HH 

891 


NL= 1 

892 


IF C IPART-1 ) 101* 10 2* 101 

893 

101 

IF C LIN ) 103* 103* 102 

894 

102 

NL = 0 

895 

103 

IF( IP*ART-1) 104* 104*15 

896 

104 

IF(LI-N) 105* 105* 15 

897 

105 

H= HH/2tE0 

898 

4 

IPT = IG0*NRYO 

89 9 


D3 5 I * 1 * NVEC T 

900 


IZ « IPT + I 

9 01 

5 

C ( I Z ) = Y(I) 

902 


igo=igo+i 

903 


I F ( IG3-2 ) 15* 15* 1 

904 

15 

03 20 I=1*NVECT 

905 


IZ = I SUB + I 

906 


Y K ( I ) s= Y d Z ) q 

YPd) « 0. 

CONTINUE ,4 

H06 » H/6 * 

Bd) » H/-2 » ^ 

B ( 2 )'=B t 1 ) 

B(3)=H 

907 


908 

20 

909 


910 


911 

912 


913 


TN( 1 ) *T I 

914 


T N 1 2 ) = T I +B ( 1 ) 

915 


T N ( 3 ) = TN ( 2 ) 

916 


TN( 4 ) = T I +H 

917 


JSUB=-1 

918 


00 30 J = 1 * 4 

919 


I F ( J-3) 106* 107* 106 

920 

106 

JSUB= JSUB+1 

921 

107 

IPT = JSl)B#NR YD + 1 

922 


CALL RHS (YK*C( IPT)*C(IPT)*TN(J) ) 

923 


IF{ 0-4)25*50*25 

924 

25 

DO 30 1 = 1* NEQ 

925 


YP( I )=YP{ I ) +A ( J )*LF ( I ) 

926 


IZ = I SUB + I 

927 

30 

YK(I) = YdZ) + B ( J ) *LF ( I ) 

928 

50 

DO 60 1 = 1* NEQ 

929 


I Z= I SUB+ I 

930 

60 

YdZ) = YdZ) + H06*{ YPd >+LF d ) ) 

931 


IF < IPART-1 165*65* 1 

932 

65 

I F t LIN) 68*68* 1 

933 

68 

T I=TN ( 4 ) 

' 934 


GO TO 4 

935 


F NO 

936 


SUBROUTINE RHS< Y* YO*Y0*TN> 

937 

C 


938 


REAL Y0(20)*YO(20)*LF*PY(20)*U(20)*Y(20) 

939 


COMMON F 1 C < 40 ) * P ( 20* £ 1 ) * C ( 80 ) * LF ( 20 ) * S ( 2 ] * 22 ) * 

940 

1 

RD1( 9) *TRACE*RD?(5)*H*RD3(2) *T*RD4 I7)*TBV( 21 )*BV(80)* 

941 

2 

AUX-< 10)* CO (67)* 

942 

3 

101 ( 16) * IPART* JSUB* 102(2 ) *NL* I03{ 29 ) *NRANKC* ID4 ( 62) * IQ( 24 ) 

943 


EQUIVALENCE ( LF * D Y ) * < P* U 1 



944 
94 5 

946 

947 

948 

949 

950 

951 

952 

953 

954 

955 

956 

957 

958 
.959 
960 

J3LLL 

962 

963 

964 

965 

966 

967 

968 

969 
97Q 

971 

972 

973 

974 

975 
975 
977 
97 8 

979 

980 

981 

982 

983 

984 

985 

986 

987 

988 

989 

990 

991 

992 

993 

994 

995 

996 

997 

998 

999 
i.000 
1P01 
1<?02 
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NRC » 20 

IF< NL > 45, 4 b> 5 0 
C 

C HERE ARE THE NGN-LINEAR OR LINEAR) EQUATIONS 

C 

45 CONTINUE 

DY< 1 ) = Y (2) 

D Y ( 2 ) * -Y( 1 ) *Y { 4 )-Y( 2 ) *Y ( 3 ) + SIN( Y<5 ) *TN ) 

47 RETURN 
50 CONTINUE 
C 

C HERE ARE THE LINEARIZED EQUATIONS 

C 

LF{ 1 ) = Y ( 2 ) 

LF(2}=-YO(l)*Y(4)-YO(2)*Y(3)-YO(3)*(Y(2)-Y0(2n-YO<4)*<Y(i)-Y3(in 
1 +SIN(Y0<5J*TN) 

RETURN 

EMI 

SUBROUTINE WRITES 
INTEGER I JT V ( 20 ) 

INTEGER QBV,XACTIC, SKIP, I B 1 G ( 135 ) , X AC TP V 

REAL RBIGI204 >, LOWER, SV.EC{ 462 ) , LF, IC, I IC, C ( 462 ) , Y ( 20, 21 ! 

COMMON IICI20),IC(20),YVEC( 420), YNU 80 ) , LF ( 20 ) , S ! 2 1 , 22 > , 

1 DELT, 5PTRB, PNflRM, TST AR T, TSTOP, 

1 OET, SMALL, ZERO, DET-T, TRACE, 

2 DZERO, CON V, DPTRB, DTEMP, DN 1 , 

3H,RN1,RN2,T, 

4TAVG,TL,TN,TBMT,TEM,TEMP,TTEM,TBV( 21 ) , BV ( 20), 

5 UPPER < 20 ) , LOWER ( 20 ), PTRR (20) , AUX { 10},CQBV< 20 ), 

6DIC (-20)-, YTPRB ( 20 ) , DEL TP, DEL TM, TNORM, S NO RM, ANORM, T BVK, DOT, 

1 IOUT, LIN, IQBVS, LITER, NORMRD, 

1 ■ NBV, NEG, NCN, NWRITE, SKIP, 

1 ITRACE, LIMIT, I VAR Y, IAF(2), 

1NVECT , I XTRA, 

1 JXTRA,KXTRA,N6UT, NL,NSPACE,NRCD, 

2NRSD,NRYD, I, II, ICOUNT, INI, IM2,INTKEY, IPLJNT, 10, IS, I SUB, I TEMP, ITER, 

3 J, JJ, JQ, JSUB,K,KK, KOUNT, KOUT, KQ, KTEM , M I TER, NB VM, NC OL , NR A NKC, 

4 NR0WS,QBV(2O),XACTBV(21),XACTIC(2O},IQBV(?O),NR,IGO,NE,NI 
EQUIVALENCE < YVEC,Y ), (S,SVEC) , ( YVEC,C), ( IOUT, JBIG), (DELT, KB IG) 
EQUIVALENCE (LF,IJTV) 

C WRITES SHOULD NOT BE ENTERED WITH LIMIT = 3 

GO TO ( 172, 200, 200, 230), LIMIT 
172 IF( IOUT-1 ) 182, 174, 174 
174 WRITE (6,550) 

DO 180 I=1,NR0WS 

180 WRITE (6,-800) I , S ( I , NCOL ) , ( S ( I , J ) , J= 1 , NR A NKC ) 

182 RETURN 

200. IF (IPUNT-2) 201,205,201 

201 CONTINUE 

IF (IOUT-1) 182,182,202 

202 DO 203 1=1, NR ANKC 

203 WRITE (6, 360) I , YNL ( I + 60 ) , I JTV ( I ) 

WRITE (6,650) ( S ( I , NCOL), 1=1, NR A NKC) 

211 WRITE ( 6, 370)D0T, ANORM 

WRITE (6, 380) ( D I C (' I ) , I =1 , N VECT ) 

RETURN 

205 LIMIT = 3 
345 WRITE (6,340) 

RETURN I 
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1903 

1904 
1005 

1906 

1907 
1008 
1909 
1019 

1911 

1912 

1913 

1914 

1915 

1916 

1917 

1918 
1019 


230 M»1 

IF (I0UT-3)232*232* 231 

231 IF (ITER-MITER) 234*232*232 
234 M = NRANKC 

232 D0 233 J=1#M 

233 WRITE (6*700) 1 AF ( I GO ) * TL* ( Y ( I * J ) * I = 1* NPSUT > 

RETURN 

340 F ORMA T ( 44 H PUNT... ILL CONDITIONED IPUNT RETURNED AS 2) 
360 FORMAT! 10X* I5*8HTH P I VO T* E 1 5 . 6* 1 6 > 

3 70 FORMAT ( 1 X/3X * 4 HD 0 T * E 1 5 • 7* 1 5 X* 7H AN0RM* El 5 • 7 ) 

380 FORMAT (/*11H DIC VECT0R*/(1X* 7E15-7)) 

5 50 FORMAT (IX// 2X*11HSAVL MATRIX*/ 

1 5X* 46HB0UNDARY VALUE* UNPTRBD S0LN* PTR8D 50L NS 

650 FORMAT ( 1X/3X/1 IHC0NST ANTS * 7E 1 5 . 7/ ( J4X* 7E1 5 • 7 ) J 
700 F0RMAT ( 1 X* A2* 8E 1 5 . 7/ I 20X* 7E 1 5 . 7 1 ) 

800 FORMAT! 3X* 12/ 7E14«6/(5X* 7E14,6)) 

END 
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EXAMPLE APPLICATION 

To illustrate in an explicit fashion the use of the program 
and, in particular, the programming of the governing equations, • 
consider the damped harmonic oscillator with an impressed forcing 
function, 

x(t) + yx(t) + £x(t) = sin wt 

the second order equation is converted to a system of two first 
order equations by the transformation 

y x (t) = x(t) 
y 2 (t) c x(t) 

The resulting system .s seen to be 

y^U) = y 2 ( tl 

y 2 (t) = -sin uit - py 2 (t) - SyjCt} 

For purposes of identifying the parameters y, E, and w from 
an observed response of the system to a given excitation, let 
us denote these constants as variables, i.e., elements of a 
state variable vector. 
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It follows that those "variables" obey the equations 

y 3 (t) * o 

y 4 (t) * 0 

y 5 U) = o 

Note that: 

a) if y, £ and to are specified then the system 

*1 * y 2 

y 2 “ sin 'y2^ t ^ “ y 4 f t 3*y I (t) 

y 3 = 0 y 3 (0) * M 

y 4 = o y 4 co) * c 

y s = 0 y 5 (0) = to 


is linear or nonlinear, as you wish to call it, 

b) If say, y and E, are not specified, then the same 

system is nonlinear for the identification problem. 

In such a case, the linearized system is constructed 

from a previous approximation y to be 

o 

• . 2 

y 2 = ^ a 2 j(t) -y^-Ct) + f 2 ct) 

j = i 

y 3 * 0 

y 4 - 0 

y -s 3 


o 



where 


4 

. f 2 (t) = sin uit - y (t) -y^ctj - y^Ct) -y^Ct) - ^ (t) • y Q (t) 

j “1 

with 

a 2 j(t) x syj £ sin ut - y- 3 ( t )’y 2 Ct) - y 4 (t) -yjC-t) 3 y * y o 

(5.9) 

such that 

a 2i(t) - -y 0 (t) 

4 

a 22 ( t ) = ~ y oJ‘ t ^ 
a 23 (t) = ' y o 2 (t) 
a 24 (t) = 

The nonlinear system can be programmed as 
D Y ( 1 ) = Y ( 2 ) 

DY(2) = SIN (Y (5) *T) - Y(2)*YC3) - Y(1)*Y(4). 

The linearized system can be programmed as 

LF ( 1) * Y (2 ) 

L F ( 2 j = SIN (YO (5) *T) - YO(2)*YO(3) - Y0(4)*Y0(1) 

+ ( Y (1 ) - YO ( 1 } ) * ( -YO (4) ) 

+ (Y(2) - YO (2) ) * ( - YO (3) ) 

+ (Y(3) - Y 0 ( 3) ) * (-YO(2)) 

+ (Y(4) - Y 0 ( 4 ) ) * (. - YO (1) ) . 
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*1 N P u T ’ Da t a "cah'o 
I N P U T_ D A T A CARO 
"INPUT DATA 'CAR"0 
INPUT DATA CARO 



1 NPUT D ATA_CAR 0 1 Q. 2 0 00 0 p O £ - 01_ 

"INPUT DATA CARD 1? b.9VS999S‘E-0A 

__ BLA'NK data CaRu 


Input oat/T'card i '" '"0 d. iuooodoe+oi"’-o,2?oooooe+oo ' l 

INPUT_DAT_A card 2 p 0.2JOOOUOE + U1 0.3S30000E-01 1_ 

INPUT DATA""cARD 3" 0 0.3UD0l’DDt'*Ul'"-0’. A 7 AOOOOE +00 1 


Input data card a _ 'o _ v ,. -.oooooot+ui -o.Ssgoc ode+oo 

INPUT DATA CARo' ' &' 0" U.bOuOOOO£+Ol' O . 3929999E+00 

INPUT Data CARD 6_ _ 0 u .bOOOOOOE+Ol 0. 1596999E+03 

In'PUT, DATA CARD""" 7" '0"‘u.7O00000£+Ul 0. 1452000E + 03 

Input data c*ro 8 o o.soocoooE+oi -o.3«8opoot+oo 

InPuT OATA"CAHD " 9 0 J.iOOCOOvE + oT -072323999E+0J 



INPUT DATA CARD 10 

InPuT oata'‘cako “if 
input DATA CARD_ 12 
Input data card' “13 
input oa_ta_card_ 14 
InPuT'data car'd lb 


0 0. 1DODU0OE+O2 -0.2273999E+03 _1 

\j “u. 1 3 000u0£+u2 0»8 hOOOOOE-01 3 

0 0 • 120U000L+02 0 . 27 1 0S99E+0 3 1 

0 u* 33OOO0OE + U2 0. 29970C0L+01 1 

0 U • 3 A 00000 E + 02 _0 « AO 1 OOOOE+OG 1 

(7 o" 1500000E + U2' -6'.28i6OOOE + 0T V 



BLANK DATA CARD 


"InPuT Data CARD ""f " 0 "d • bOOOOOOE + UU 0. OOOOOOOE+OO "“CVnonocOOE+OT; '0. 00060002 +‘0t)‘ 

INPUT Data CARO S 0 >, . iOOQCOGE + 0 L O.OOOOOOOE+OP 0 • OOOOOOQ E + 00 0 • 00 000 COE + 00 

"InPuT DATA CARD 3 0* >/• 99999 9oE- 01 0 .OOUOOOOE+OO 0“'f)000'006'E +00 - 0 . Od'OOOOOE + <J0" 

INPUT DATA CARD A 0 J. .299999E+01 O.OOOOOOOE+OP, 0. OOOOOOOE + OO 0 . OOOOOOOE+OO 

INPUT Data CARD b -1 o. lOOOOOOt + Ol 0. OOOOOOOE+OO 0 • OOOOOOOE + OO "O.OOOOOOOfc+OO 
BLANK DATA CARD ... „ .. 


TBV ' B-V QBV XCT8V 

1 0 . lOOOOOOE+Ol "0 • 22DOOUOE + 00 0 1 

'2 "0.20000UDE+01 ' 0.3500UUPE-01 0 ‘ 1 

3 0. 300060UE+01 -O-.A/aOOOOE + OO 0 1_ ^ 

“a “ d.AOOC-O'OOt + Ol -0« bR90vU!)F +00 “ "o’ f 

5 O.bOCOOOCE + Ol _ 0*39299-j9L+00 0 _ 1 

‘ 6 ~" ' 0 • fcU000uUL+01 0. 1596999F + 01 " ' 0 "" ""l 

7 O.VOOCOOCE+Ol 0 . 1 Ab20oCF + 0 1 0 1 

"a ~0. SuOOOOOE+Ol' “0 • 3880 OOOE + 0 0 O' l‘ 

_ 9 0 •_ 9 u 0 0 0 0 0 E + C1 -0_. 2J23999E+01 0 l 

10 0 ."ToOOOOOE + OE “0 * 22 7 J399E+01 'o 1 

11 0 . 11 OOOOOE+02 0-8300000E-01 0 1 

12" 0. 12OCOD0E+-O2'" “0.27109992 + 01' " 0 T 

13 0 . 1 300000E+02 0 • 2997U00E+0 1 0 1 

14 6". 1 7000002 + 02 ~0* AOIOOOOE + OO' *6 1 

15 .0, lbOOOODt+02 -0 • 281 6000E+01 0. 1 



n vnb & 




y< ) 

•1C 

XCT1C 

PTRB 

UPPER 

LtiwER 

1 

o.sooooodi+oo 

o' 

0* JPOOOOOL+Ol 

o.ooOooocE+oo 

O.OGOO000t+00 

2 

0. lOOOOOOE+Ol 

0 

6* lPOOOOOE + Ol 

O.OOOOOCOE+OO 

O.OCOOOOOE+OO 

3 

0. 9999996 E“Ol 

' 0 

0* 1200000E + 01 

o'. o'ooooo'Oe.+Oo 

O.O.CO’uOTBeT+TH) 

4 

0*1 299999E+01 

0 

0* 1 200000E+01 . 

O.OOOOOOOE+OO ■ 

0 < OOOOOOOE+OO 

5 

0* 1000000E+bl 


'0 ■* 0 0 00000 e "+ 00 

OVOOWCOOE + OO 

"0.5000()0"OE + 00 


• ieur o ' ijn DI Q3VS o liter go norrrp, o ; ■ 

N'BV 1 5 ME a 2 MCN 3 NWRIIN^ j! SKJj=> ii . 

P ELT O.gCOOQOE-01 SPTRB 0 . 1 20000E+0 1_ PN0RN 0._OOOOOOE+O0 TSTAR T Q.O OOqOOE+OO TSTtlP 0 » IbOQQOE+Og 

SAMPLE GUT POT ~ . 

GENERATED 

NRtfNKC 5 . 

AUX " ' ' " 

o. ooocoooe+oo _o*oooooooL+<>o d* oooouooe+’oo o*ooooo do e+oo o.oooooooe+oo o * oooooooe+oo ooooocof +oo o*oooooooe+oo 

'"'0"*0006000E + 00“ "0*00000000 + 00 ‘ - ■ “ 

~ ITER I ilr 20 

..TIME ■*“ ~y'OLUfTfiN 

0<OCOOOO OE* QO Q.SOOOOOOE + OO 0 . lOOuOOOE+O 1 0.9999996E-01 0 * 1 299999E + 0 1 _ 0 *J OOOOOOE + Ol 

‘ ITER ' 2 l>F 20 ' ' ' ‘ ' " ; 1 

TIME SOLUTION 

pj'doooocof+oo ' 6*834114 1 E. + 00 -0*S40l550E-01 6 .*36 7'2? 28 E + OO 0 i'i 158999E+01 “ 0* 1 OOtSOO'OE+O j‘ ‘ ‘ ----- — - 

ijLR__i jf. 20 ; .. 1 

tine ’ Solution 

P*OOOOOCOE+O0 -0*£206699t. + 00 0 . 2298086E + 01 0 . 1 94 9387E+' 00 0.86186182+00 0 . 1 00 OOOOE + 0 1 

‘ ITER 4 tiF 20 ' ‘ " - T 

-TINE SOLUTION 

. CJOOOOOOOE+OO -0 »5682ab9t+00 0 . 32b6903E+00 072'36625'7E+00 Ti~9734 84 8E + 00' " 0 .TOu'CTOOOE+Ol 

.. ITEK_ b .OF 20 . __ _? 

"fl>E • " SOLUTION 

O.OOCOOOOE + OO 0 • 1 36 25 + 6E + 0 1 U . 3 < 76879E + 00 0 . 201 7888E+O0 0 . 1000775E + 0 1 0 . 1 OOOOOOE + Ol 

. - .. mR h ,. F ?0 ~ ' " ' ■» 

t i ^ e solution _ _ _ _ _ 

O'SoOOObOEJuO C *998b23 7E + 00 0.'5017?81E+6o“ oVl 9990 71 E + 00 6'.“995'8'708E + 00 'O'* 1 OGCfoOOE + 01 ' 

ITER 7 t’F 20 . j 

TIRE • "sOLUT I BN ” 

0 » OQ OOOOOE + 0 0 0* 100202, ,-E + Ol 0 . 50 21 760E + 00 6 « 1 999351 E+00 0 •9998615E + 00 0 » 1 OOOOOOE+Ol _ __ 

ITER “ 8' "Jr ~20 ” " " ““ T 

TIRE SOLUTION ' 

O.OOOOOCOE'+OO 0 . 10019 />t+‘Ol" 0*'5O2U00 Ve + 0O 0“« i 999393E+00 0.99986S8E+00 0 * 1 0000'00£ + 0'l 

■ . ^ I 


u 

■p 



Iter V '"'f 9 

time _ _ ; S0LUTIBN _ j 

0«OOdOOOdE+UO oVrobiaEj'+Ol” 0.50?br53F+00 o 

Q. 99999965-01 0. i 04 6 7506 + 01 0 . 39552 73 E + 00____ 

0» 1999999L+"0O 0*108134/ £+01 “ 0 . 89 7o799E +00 

0*29999y8E+00 0 . 1 1 06 4 o46 + 0 1 C * 207 j300E + 00 

0* 3999597E+00 0 . 1 l23l'l l E + 6l“ 'o . 1 26 748 7E + 00 

Oj 4999996E+00 0_* ll32l4bE+qi 0.55c 1259F-01 

0 • 5999995E+00 0 -’ll 3454>6 + 01 -0 .599 1302E-02 

0 * 6993994E + 00 0 • 1 1 3 1 251 E + 0 1 -0 . 58,1 591 OE-0 1 - 

0. 7999393E+00 'O • 1 1 23?0bt+0l -0. 1 01 1595F+00 

_ 0.S999992E+00 0 . 1 1 1 1 JU'E +01 -0 * 13542?3E + 0C 

BV 6. 1000Ga'ut + 01~'o'. 10963836 + 01 -0 . 1 6 i b?4 3E+0 0 

' 1TH BV >0. 2 80000 024-00 i, (Jr ) -0 . 2224 722E+00 

0‘7 1 0 9 9 9 9 8 E‘+ 6 1 O'V 1 b 7 9 2 34 E + (> I" - 0 . 1 8 U 1 7 2 2 E -i- 00 
O'* 1 1 99995E+Q1 0 . 1 06056<-£ + 0l • -0 . 1 92 1865E + 00 

0+ 1S99993E 401 6. 104 0 979t4 01 -0 . 1 98 + 8 l 7F + 00 ' 

0* 1399930E+01 0. 102101 3 E+01 - 0 . 20 00 46 1 E 4 0 0 _ 

0. l«99988E+01 ' 0. IOOIOS'-E + Oj ‘ -0.19/3221F+00 ‘ 

• 0. 15999S5L+01 0*961 5i i 5E+00 -0 . 1 931 84 5>- + n0 

'0^16999832 + 01 ‘O- 96249 7cE4JO“'-0.18e9200E+o'0 

0. 17S9SR0E+01 0.94414192+00 -0 . 1 HO2057F+00 

0- 1899978t40r'""0. 9264361 E400 - 0 . 1 7 4 0 88 8E + 0 0 
BV 0« 2000CGOE+01 0 . 9092662E+00 -0 . 1 69b670F+00 

2TH flv 0.3500000E-01 OIY > 0 . 340541 OE-01 " 

_____ 0 • 2C999 1 .! 82 + 0 1 _0.89243dl!t400 -0 • 1 6 75733E+00 

“0*213999bE4dr‘ 0". 8 7563 8nE + 00 '-0 . 1 689528E + 00 

0 » 2299593E -r 0 1 0 * 8=«50cl"i + 00 -0 . 1 744520E400 

“ 0. 2359999E + 01 " 0 * 84Ub85( 2 + 00 ' -0 . 1 84 703oE + 00 

0.2699988E + 01 0 . 821 38*4 fc. + 00 »0 . 2002 1 4 3E + 0 0 

6. 25.39985E + 01 “ 0 . 8003b /5E+00 -0 * 221 J57 DE + OQ' 
0.26999H3E+01 0 • 7 7 6920 /E+00 -0 . 24 236? 1 E+00 

0 * 2799 98OE + 0 i 6 . 7 5048b?fc + {7C -0 . 2813 123E + 00 

0 . 28999 7SE + 01 ’ 0 • 7204 6 1 6E +00 -0 . 3201 4 1 1 E + OO 

“BV “O- 3009000E+01 O'* 6 8 o2bi»4fc + 00 -0 * 3 6 4'o4 3 3£ + 00' 

; 3TH BV -0. 4 /400C0E+00 u(Y) -0 . 4/2 1 429E +00 __ 

' 0 . 3C99998E + 0 f ‘ 0 . 6 4 734 7 It + 00 -<*0 • 4 1 4 4333F + 00 

0. 31 99 9 3 5E +01 0 . 60321 2t E + 00 -0.4 6 901 3 9F + 00 

0 •" 32999932 +01 (f* 553 406 1 b + 00 -0.b2 7/4 1 9E + 6cf 

0. 3399990b + 01 0 • 4 9755 1 t + 00 -0.589 84 7 46 + 00 

0.34995882+01" 0 • 4 3b3530t + 00 -0 . 654 444 6E + 00'“ 

0 * 3599985E + 0 1 0.36661 14 6 + 00 -0. 7205452E+00 

0.3639S83E + 01 ' ‘ 0 *29 1 229CE + 00 -0 .78 /0732E + 00 

0. 3739 9808 t01_ 0*2092201 E + 00 -0 * 85288 1 OF + OO 

' 0 * 36999782+01 o’« l 207 1 64 2 + OO" - 6 . 9 1 6 76 7 2F. + O0 

BV 0. 400COQOE+01 0*25347892-01 -0 . 977509 1 F+00 

4TH Bv -0*58900002+00 l, i Y > »0 ■ 58 730492+00 

0 * 4.0999S8E +01 -P • 7466060E-0 1 -0 * j 03382 1E + 01 

61+1 S99S5E+G 1 -0*180627. 5 6 + 00 -0 . 1 0844 72F + 0 1~ 

0 * 429993'3£ +01 -0 . 29 1.324 /t +00_-0 . 1 1 26238E + 0 1 

0. 4 39399C.fi ♦Of -0. '406005 46 + 00 -0*11639381+01 

0*44939882+01 -0* 52380566+00 -0 . 1 1 90456E+01 

"0.4593985E+0i"-d.6437bb2fc+Od -0 * 1206755E + 0 l r 

. 0 * 4699383 E + 01 - 0 . 764 784 22 + 00 -0_. 1 2 US02E + <M 

0.47999802+01 -0 • 8ab 7359E + 00 -0 « 1 2U5082E+0 1 

0. 489 99/B E+Q1 -0 ♦ 1 OOb37/iE + 0 1 -0.1 1656 10E+01 

BV 0.50000 00 £+01 -0.11224426+01 -0 . 1 13294 5E+01 
t>TH BV 0 • 3929999E +00 0(Y ) 0 * 393S88 4F + 00 



(. 


_ _p.50939^88+pl_ -C..1235535E+01 -0 . a 1 06721 c +01 _ 

’ Oi5l9999b£ + Ol -0* 1 3,433188 + 01 ~ 0 . 1 0*6728E+oi 

p« 5299 99 3t+Q 1_ -p • 1 4 4 4 + 12£_+ 0 1_- 0 .97294 1 5£ + 0_0_ 

0.5897990L + 01 -0 . 1 537V*<t 8 + 01 -0 . 8855 1 30F + 00 ' 

, 0.5*999831+01 “0 . lb2106*E+01 -0. 784 78?*E+00 

0.559998bE+0i‘-d. l693‘ybK,E + 01- -0.6712782E+00 ' 

0.5699983E+01 -0.17549072*01 -0 . 54 5/ 1 52E+00 

• 0 • 5739950E +0 i" -0 . 1 80272 /t + OT -O'. 4089893E + 00“ 

P+5899978E+01__^p. 18363632+01 -0.562171 7F+00_ 

BV 0.5q 60O0OE’+'01 -0'.1854 862E+01 -oil06461?E+'nO 
6TH BV 0. 1596999E + 01 Ci(.Y) 0 * 159848 7E + 0 1 
"0 . 60999S8E+01 i;0« 1857+032+01 0 • 56e7 4 b3E-0i~ 

_ 0.6199995E+01 -0 . 1 843321 2 + 01 0 .2,257376£ + 00_ 

0.629 9,9 93E+bi"«0»i8l210ot + 0i ‘ 0*39907 1 3E+0Q 
0.63999902 + 01 -0 . 1763* 1 7E+ 0 1 0 • 57490s’7E + 00 

0'.6*93388E + dl -Q.' 1697096t + 0l 0. 7bli9*8E + 00 
__ 0.65999S5E+O1 -0 . 1 o 23 1 73E+01 ’ 0 . 926o07 1 F +00 
0.6699983E+01 ~-C» 151187 jfc+01 '0.109357*8+01 
0 . 6799980E+01 -0. l3«36238+0l 0 • 1 2oo306F +01 

, ‘ 0.6699978L+.01 -0 . 1 25904&E+ 0 1 0.142*8158 + 01" 

BV O'. 70000CQE + 01 -0.U0o92l£ + 01 _ 0 . 1 b 7a 1 69E +01_ 
77 H faV 0. 14520008+01 ' i/m ' 0. 145032 + E + Ol 

, r 0.70999988+01 -0. 94 4340 o £+00 , 0 . 1 7 1 43 79E +01 

" 0.71 9999bE+0l'‘-0. 7664 738E + 00' ,0. 184064 3F + 01 
___ 0 . 7299 993E + 01__”0 .5767014E + 00 0 . 1 <Ja2 2 1 0E + 01 

' 0.739999OE + 0l"-O.3765>b2t + 00 ' 0 . 2047459F +01 
_ 0.74999,r8E + yl_-0.16/8u2yt+Oe_ 0 . 2 1 2+ 90* F +01_ 
’ O'. 7b99935E + 0'l 0 . 4 7 7 o 3 1 1 8 “ 0 1 "C . 2 1 8 J2??E +01 

0. 7699983E+01_ C . 268 1 bbftt+ 00 0.2221?7*E+P1 
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p._2499988£+0t 0.8218322E+00 _-0 . 1 9532 1 lE+OO 

BV 0 • 3000000E +0 1 0. 6888731 E + 00 -6 . 36 1 0 1 67E+00 

Q.3499988E+01 _ 0 • 4 393052E+00 -0. 6527909E+00. 

BV 0.4000000E+01 0 . 3018703E-01 -0. 9779964E+00 

0. 4 4 99 988E + 01_«0 # 5202802 E + 00 -Q . 1 1 92725E +01 

BV 0.5000000E+01 "0 • 1 1 20350E+01 -0 . 1 1 56263E +0 1 

P * 5^ 9998 8 E + 01 *0 ».l 620703£+01___“0 • 78 BH4S8E +.00 ; 

0. 599997SE+01 »0 . 1 856094E+0 1 -0 . 1 092827E+00 

P*6499963E+0 l_ltP.* 1 699442E+01 . . 0 . 74 98 1 4 7E +00 

0* 6999950E+01 -0 . 1 1 1 173* E + 0 1 0. 1574995E+01 

. p.7499938E+Ql_-0.l703653E+00 0 . 2 1 26025E+0 1 

BV £)• 8000000E + 01 0 .93556528 + 00 0 . 220 74 1 7E+0 1 

p. 8499988E+01_ 0_. 1 943647E + 01 0 . 1 7 3256 3E +0 1 

0* *89999758+01 0 . 2584495E+0 1 0 . 76 1 1 488E + 00 

Q.949996 3E+01_ 0.26553198 + 01 -0 . 5033699E+00 

0. 9999S50E +0 1 0.20826358 + 01 -0 . 1 767638E +0 1 

0. 1049994E+02. .0.9543239E+00 »0 • 26 73S94E+0 1 

BV p. 11000C0E+02 "0 . 4902a 1 6E+00 *0 . 29K7592E+0 1 

0. 1 149999E+0 2 -0 .191 1094E+01 -0 . 25 73289E+0 1 

0.11999988 + 02 -0 .2950648E+01 -0. 1487777E+01 

p_. 1249996E + 02_.-Q.;3325378E + 01 0 . 34 06220E-0 1 

0. 1299995E + 02 -0 . 2903840E + 0 1 0 . 163078SE+01 

p._1349994E+02_-0_..l749135E + 01 0 . 29C21S6E+01 

Q. 1399993E+02 -0 . 1128750E+00 0. 351 1342E+01 

0.14 49 9S IE + 0 2 0 ■ 1 620046 E+01 0 .3 2. 7 4 6 6 3 E + OJ . 

BV p.l5po6odE+d2 0.3023725E+01 0 . 22 15425E+01 


